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(Manuscript  received  18  October  2013,  in  final  form  10  June  2014) 

ABSTRACT 

Pressure  differences  across  topography  generate  a  form  drag  that  opposes  the  flow  in  the  water  column,  and 
viscous  and  pressure  forces  acting  on  roughness  elements  of  the  topographic  surface  generate  a  frictional  drag 
on  the  bottom.  Form  drag  and  bottom  roughness  lengths  were  estimated  over  the  East  Flower  Garden  Bank 
(EFGB)  in  the  Gulf  of  Mexico  by  combining  an  array  of  bottom  pressure  measurements  and  profiles  of 
velocity  and  turbulent  kinetic  dissipation  rates.  The  EFGB  is  a  coral  bank  about  6  km  wide  and  10  km  long 
located  at  the  shelf  edge  that  rises  from  100-m  water  depth  to  about  18  m  below  the  sea  surface.  The  average 
frictional  drag  coefficient  over  the  entire  bank  was  estimated  as  0.006  using  roughness  lengths  that  ranged 
from  0.001  cm  for  relatively  smooth  portions  of  the  bank  to  1-10  cm  for  very  rough  portions  over  the  corals. 

The  measured  form  drag  over  the  bank  showed  multiple  time-scale  variability.  Diurnal  tides  and  low- 
frequency  motions  with  periods  ranging  from  4  to  17  days  generated  form  drags  of  about  2000 Nm“^  with 
average  drag  coefficients  ranging  between  0.03  and  0.22,  which  are  a  factor  of  5-35  times  larger  than  the 
average  frictional  drag  coefficient.  Both  linear  wave  and  quadratic  drag  laws  have  similarities  with  the  ob¬ 
served  form  drag.  The  form  drag  is  an  important  flow  retardation  mechanism  even  in  the  presence  of  the  large 
frictional  drag  associated  with  coral  reefs  and  requires  parameterization. 


1.  Introduction 

Roughness  elements  and  rough  bathymetric  features 
generate  resistance  to  flow  through  skin  friction  and 
pressure  anomalies,  which  are  referred  to  as  frictional 
drag  (or  frictional  skin  drag)  and  pressure  drag  (or  form 
drag),  respectively.  Frictional  drag  is  a  tangential  stress 
on  the  surface  caused  by  the  molecular  diffusion  of 
momentum  across  a  velocity  interface  (Schlichting  1962) 
and  also  applies  to  heat,  salt,  and  other  scalars.  This  is 
primarily  true  for  hydrodynamically  smooth  flows  where 
the  bulk  of  the  force  on  the  surface  is  by  viscosity  be¬ 
cause  pressure  forces  are  negligible.  However  in  rough 
flows,  where  the  roughness  elements  are  larger  than  the 
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molecular  sublayer,  the  force  on  the  surface  is  transmitted 
mainly  by  the  pressure  forces  on  roughness  elements,  and 
the  force  is  parameterized  as  frictional  drag  (Schlichting 
1962;  Tennekes  and  Lumley  1972).  The  frictional  drag 
depends  strongly  on  the  hydrodynamic  roughness  or  the 
roughness  Reynolds  number  (e.g.,  Nikuradse  1950). 
The  parameterization  of  frictional  drag  is  typically 
derived  from  inertial  or  wall  boundary  layer  theory  (e.g., 
Tennekes  and  Lumley  1972;  Stull  1988)  and  has  been 
applied  in  both  atmosphere  and  ocean  models  (Blumberg 
and  Mellor  1987;  Hodur  1997). 

Form  drag  is  related  to  the  dynamic  pressure  differences 
across  an  obstacle  arising  from  upstream  blocking  and 
separation  of  flow  and  formation  of  eddies,  hydraulic 
jumps,  and  lee-wave  disturbances  downstream  of  the  ob¬ 
stacle  (Gill  1982;  Smith  1989;  Baines  1995).  The  nature  of 
the  flow  depends  on  the  governing  parameters  such  as 
Froude  number  Fr,  nondimensional  obstacle  height  hn  and 
width  dn,  and  mode  number  K.  For  single-layer  hydraulic 
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flows,  Fr  =  where  Uq  is  the  flow  speed  of  the  layer, 

g  is  the  gravitational  acceleration,  and  D  is  the  flow  depth. 

In  stratifled  depth-limited  flows,  the  nondimensional 
height  of  the  bank  or  nonlinearity  parameter  /z„,  the 
nondimensional  width  the  mode  number  K,  the  aspect 
ratio  of  the  bank  ho/L,  and  the  height  parameter  ho/D  are 
used  as  nondimensional  parameters  to  describe  the  na¬ 
ture  of  the  flow,  where  hn  =  NoholUo,  an  =  and 

K  =  NqDI{ttUq)\  No  and  Uq  are  characteristic  scales  of 
stratiflcation  and  flow  speed  of  the  water  column,  ho  is  the 
height  of  the  bank,  D  is  the  flow  depth,  L  is  the  width  of 
the  bank,  and  ao  is  the  half-width  scale  (Hunt  and  Snyder 
1980;  Smith  1980;  Durran  1986;  Baines  1995;  Vosper  et  al. 
1999).  For  hn<l  and  K<1,  the  flow  moves  over  the 
bank  (nonblocked  flow),  and  for  hn>l  and  K>1,  the 
flow  moves  around  the  bank  (blocked  flow)  (e.g..  Hunt 
and  Snyder  1980).  In  the  presence  of  background  strati¬ 
flcation,  the  flow  over  obstacles  generates  internal  gravity 
waves  and  in  turn  transports  horizontal  momentum  ver¬ 
tically  (when  hn<l).  The  resulting  wave  momentum  flux 
generates  a  form  drag  to  the  background  flow  that  is  also 
referred  to  as  wave  drag  (e.g.,  Durran  1986, 1990;  Smith 
1989).  Form  drag  is  patchy  in  space  due  to  the  spatial 
variability  of  the  bottom  topography,  unlike  frictional 
drag  that  occurs  on  the  entire  bottom  surface. 

Form  drag  for  airflow  over  mountains  has  been  stud¬ 
ied  extensively  in  both  held  observations  and  numerical 
simulations  with  special  emphasis  on  improving  the  pa¬ 
rameterization  of  mountain  wave  drag  (Smith  1978, 1979, 
1980;  Bougeault  et  al.  1990;  Durran  1990;  Clark  and 
Miller  1991;  Olafsson  and  Bougeault  1996;  Doyle  and 
Jiang  2006).  Several  decades  of  mountain  wave  studies 
were  beneficial  to  the  estimation  of  mountain  wave  drag 
based  mainly  on  linear  internal  wave  solutions.  Such 
parameterizations  have  resulted  in  substantial  improve¬ 
ments  in  current  weather  prediction  models  (e.g.,  Kim 
and  Arakawa  1995;  Wood  et  al.  2001;  Kim  et  al.  2003). 
Form  drag  has  been  computed  for  both  small  and  large 
atmospheric  mountain  ranges.  On  a  smaller  scale.  Smith 
(1978)  evaluated  pressure  drag  in  the  Blue  Ridge 
Mountains  in  the  United  States  using  several  micro¬ 
barographs  across  the  ridge,  resolving  the  pressure  field 
over  the  ridge  similar  to  the  oceanic  observations  made 
by  Warner  et  al.  (2013).  However,  for  most  atmospheric 
observations,  measurements  of  pressure  over  mountain 
terrains  are  quite  sparse,  as  the  instruments  cannot  cover 
the  entire  mountain  range.  Hafner  and  Smith  (1985)  used 
a  limited  number  of  microbarographs  over  the  European 
Alps  to  evaluate  the  pressure  drag  by  assuming  lateral 
pressure  gradients  were  constant;  even  though  these  ob¬ 
servations  do  not  resolve  the  entire  pressure  field  over  the 
complex  terrains,  a  rough  estimate  of  form  drag  was  ob¬ 
tained.  In  recent  mountain  wave  experiments  such  as  the 


Terrain- Induced  Rotor  Experiment  (T-REX;  Grubisic 
et  al.  2008)  microbarographs  measurements  were  used, 
but  a  few  microbarographs  measurements  are  not  suffi¬ 
cient  to  compute  the  pressure  drag  directly.  Instead, 
pressure  measurements  have  been  used  to  obtain  in¬ 
formation  on  internal  lee  waves. 

In  the  ocean,  observational  studies  of  form  drag  are 
very  limited  (e.g.,  Nash  and  Mourn  2001;  Warner  et  al. 
2013).  Form  drag  acts  on  the  entire  water  column  over  the 
topographic  feature  and  leads  to  the  generation  of  mixing 
and  turbulence,  eddies,  and  internal  waves  away  from  the 
bottom  (e.g.,  Polzin  et  al.  1997;  Jayne  and  St.  Laurent 
2001;  Pawlak  et  al.  2003;  Garrett  and  Kunze  2007).  Some 
aspects  of  topographic-tidal  flow  interaction  and  for¬ 
mation  of  lee  waves  have  been  studied  analytically  and 
numerically  using  nonhydrostatic  models  (Bell  1975; 
Nakamura  et  al.  2000;  Khatiwala  2003;  Legg  and  Huijts 
2006).  Several  numerical  studies  have  been  conducted  to 
examine  the  dynamics  of  flows  over  small  banks  (e.g.. 
Lamb  1994;  Skyllingstad  and  Wijesekera  2004;  Seim  et  al. 
2012).  Lamb  (1994)  and  Skyllingstad  and  Wijesekera 
(2004)  reported  qualitatively  comparable  results  to  the 
theoretical  and  laboratory  results  of  Long  (1955)  and 
Baines  (1979)  for  a  similar  dynamical  parameter  range  (1 
<  K  <  2  and  <  10).  Skyllingstad  and  Wijesekera 
(2004)  noted  for  K  <  1,  a  substantial  reduction  in  the 
wave  drag  compared  to  a  no-slip  boundary  when  bottom 
roughness  was  introduced.  Ocean  circulation  models  do 
not  have  form  drag  parameterizations  (Blumberg  and 
Mellor  1987;  Oke  et  al.  2002),  mainly  due  to  lack  of  in¬ 
formation  on  the  magnitude  and  spatial  distribution  of 
form  drag,  although  some  attempts  have  been  made  to 
approximate  wave  drag  from  linear  theory  (Nikurashin 
and  Ferrari  2011). 

The  total  form  drag  consists  of  different  components 
that  involve  isopycnal  displacement  (referred  as 
“internal”)  and  surface  tilt  (referred  as  “external”;  e.g., 
McCabe  et  al.  2006).  Different  components  of  form  drag 
over  coastal  banks  and  ridges  were  examined  by  several 
investigators  (Mourn  and  Nash  2000;  Nash  and  Mourn 
2001;  Edwards  et  al.  2004;  McCabe  et  al.  2006).  Warner 
et  al.  (2013)  were  first  to  measure  the  total  form  drag 
from  high-precision  pressure  sensors  (referred  to  as 
Ppods;  Mourn  and  Nash  2008)  deployed  over  Three 
Tree  Point  (TTP),  a  headland  in  Puget  Sound,  Wash¬ 
ington.  Warner  et  al.  (2013)  reported  that  the  form  drag 
coefficient  is  at  least  one  order  of  magnitude  larger  than 
the  frictional  drag  coefficient  at  the  bottom. 

In  the  following,  we  compute  both  form  drag  and 
frictional  drag  from  in  situ  measurements  collected  over 
the  East  Flower  Garden  Bank  (EFGB)  on  the  Texas- 
Louisiana  shelf  (Fig.  1)  as  part  of  projects  sponsored  by 
the  Naval  Research  Laboratory  [NRL;  Mixing  Over 
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Fig.  1.  (left)  Habitat  characterization  map  of  the  EFGB,  along  with  mooring  locations.  The  edge  of  the  habitat 
map  is  the  boundary  of  the  marine  sanctuary.  Habitat  characterization  was  provided  by  NOAA  (available  at  http:// 
www.ncddc.noaa.gov/website/google_maps/FGB/mapsFGB.htm).  Red,  blue,  light  green,  green,  dark  blue,  and 
yellow  regions  denote  coral  reefs,  coralline  algal  reef,  algal  nodules,  coral  community  zone,  and  soft  bottom,  re¬ 
spectively.  Barny  moorings  (M1-M8,  M4a,  and  M5a)  are  also  marked  on  the  map.  (right)  Roughness  elements 
associated  with  corals.  A  diver  was  examining  the  underside  of  a  dislodged  coral  head  after  Hurricane  Ike  in 
September  2008.  Photographer:  E.  Hickerson,  Flower  Garden  Banks  National  Marine  Sanctuary.  This  image  is  part 
of  a  larger  collection  of  sanctuary  media  found  online  (at  https://marinelife.noaa.gov/media_hb). 


Rough  Topography  (MORT)]  and  by  the  Bureau  of 
Ocean  Energy  Management  [BOEM;  Currents  over 
Banks  (COB)].  As  described  below,  the  flow  over  the 
EFGB  consists  of  /z„  >  1,  >  1,  and  Un  ~  0(100)  along 

with  large  bottom  roughness  lengths  (zo  ~  1-10  cm).  This 
parameter  regime  belongs  to  nonlinear,  hydrostatic, 
blocking  flows,  and  there  is  no  theoretical  prediction  of 
the  form  drag.  Main  objectives  are  to  quantify  the  mag¬ 
nitude  and  temporal  variability  of  the  form  drag  through 
bottom  pressure  measurements  collected  from  June  to 
December  2011  and  to  quantify  roughness  lengths  and 
associated  frictional  drag  coefficients  from  microstruc¬ 
ture  and  velocity  profiles  collected  during  an  intensive 
observational  period  in  June  2011  and  in  August  2011. 
The  paper  is  organized  as  follows:  Section  2  describes 
field  observations.  Section  3  describes  form  and  frictional 
drag  formulations.  Estimates  of  form  drags  and  rough¬ 
ness  length  scales  are  provided  in  sections  4  and  5,  re¬ 
spectively.  Discussion  is  given  in  section  6.  A  summary 
with  major  findings  is  in  section  7.  Last,  a  sensitivity  study 
of  the  estimated  form  drag  based  on  numerical  simula¬ 
tions  is  provided  in  the  appendix. 


2.  Observations 

Sampling  was  conducted  at  the  EFGB,  located  ap¬ 
proximately  190  km  southeast  of  Galveston,  Texas,  in 
the  northwestern  Gulf  of  Mexico  (GoM).  The  EFGB  is 
part  of  the  National  Marine  Sanctuary  System  and  is  one 
of  many  banks  formed  by  salt  domes  on  the  continental 
shelf.  The  EFGB  is  a  biologically  diverse  coral  reef 
consisting  of  brightly  colored  sponges,  plants,  and  other 
marine  life.  The  habitat  map  and  an  example  of  the  size 
of  a  typical  coral  head  are  shown  in  Fig.  1.  The  EFGB  is 
a  pear-shaped  bathymetric  feature  rising  from  below 
100-m  water  depth  to  18  m  below  the  surface  (Figs.  1, 2). 
The  bank  is  about  10  km  long  and  about  6  km  wide,  and 
its  projected  area  to  the  horizontal  plane  encompassed 
by  the  100-m  bathymetry  contour  is  approximately  40  km^. 
The  EFGB  has  steep  lateral  slopes  (~50-m  height  change 
within  500-m  distance)  at  the  eastern  and  southern 
boundaries.  The  top,  located  just  off  the  southern  peak,  is 
nearly  flat  with  more  gentle  slopes  westward  and  north¬ 
ward,  but  contains  rough  bathymetric  elements  such  as 
corals  (Figs.  1,  2).  The  top  of  the  bank  is  dominated  by 
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Fig.  2.  (a)  Bathymetric  map  of  the  bank  with  contours  every  10  m.  The  red  bullets  are  the  long-term  Barny 
mooring  locations,  and  the  blue  bullets  are  the  short-term  mooring  locations.  The  locations  of  Ppods  are  marked  in 
green  circles.  Note  that  the  string  moorings  (S1-S5  and  S8)  were  located  about  100-200  m  from  the  Barny  locations. 
The  green  lines  are  June  transects,  and  blue  lines  are  August  transects.  The  black  dashed  lines  shows  the  coordinate 
system  (x,  y),  where  the  x  axis  is  parallel  to  the  M2-M5  mooring  line,  and  the  y  axis  is  perpendicular  to  the  mooring 
line.  Note  that  U  is  positive  in  the  x  direction,  and  V  is  positive  in  the  y  direction,  (b)  Cross  section  of  bathymetry 
along  M2-M5-M4  line.  Distance  (positive  eastward)  is  from  M2  and  the  height  of  bank  (positive  upward)  is  from 
106-m  bathymetry  contour,  (c)  Longitude-depth  cross  sections  of  bathymetry  along  T1  and  T2  transects. 


relatively  large,  head-forming  corals,  with  Montastrea 
annularis,  known  as  boulder  star  corals,  composing  30% 
of  the  total  coverage  between  18-  and  36-m  depth,  and  the 
species  Diploria  strigosa,  Colpophyllia  natans,  Montastrea 
cavernosa,  Millepora  alcicornis,  and  Porites  astreoides 
account  for  an  additional  20%  coverage  (Fig.  1).  The  total 
coral  coverage  on  the  upper  reef  averages  approximately 
55%  (Bright  et  al.  1984).  Between  30-  and  40-m  depth, 
algal  ridges  are  common,  and  below  50-m,  algal  nodules 
(rhodoliths)  are  dominant.  In  these  deeper  zones  nodules 
can  reach  1-2-m  height  due  to  lower  current  velocities 
(Minnery  1990;  Fig.  1). 

The  observational  program  over  the  EFGB  consisted 
of  two  6-month  mooring  deployments  (from  December 

2010  to  June  2011  and  from  June  2011  to  December 
2011),  a  2- week  intensive  observational  period  (lOP) 
from  31  May  to  14  June  2011,  and  a  week-long  micro¬ 
structure  survey  during  a  coral  spawning  event  in  August 

2011  as  a  part  of  the  MORT  and  COB  projects.  Two  ships 
[Research  Vessel  (R/V)  Pelican  and  R/V  Manta]  took 
measurements  over  the  bank  during  the  lOP.  Micro¬ 
structure  surveys  were  conducted  from  the  RA^  Manta, 
while  moorings  were  deployed  from  the  R/V  Pelican. 
Detailed  discussions  of  instrumentation,  data  collection, 
sampling  methods,  and  data  processing  for  the  long-term 
and  short-term  deployments  are  in  Teague  et  al.  (2013) 


and  in  Jarosz  et  al.  (2014,  manuscript  submitted  to 
/.  Geophys.  Res.).  Here  we  describe  the  subset  of  the  data 
relevant  to  this  study,  that  is,  measurements  made  during 
the  lOP  and  in  August  2011  as  well  as  from  the  moorings 
deployed  during  the  second  half  of  the  COB  experiment 
(June-December  2011).  Bathymetry  and  locations  of 
these  short-  and  long-term  moorings  at  the  EFGB  are 
shown  in  Figs,  la  and  2a. 

a.  Moorings:  Currents,  pressure,  and  hydrography 

A  total  of  10  Teledyne  RD  Instruments’  (RDI) 
acoustic  Doppler  current  profilers  (ADCPs)  were  de¬ 
ployed  in  bottom-mounted,  trawl-resistant  housings  re¬ 
ferred  to  as  “Barnys”  (Perkins  et  al.  2000)  at  stations 
M1-M8,  measuring  east,  north,  and  vertical  velocities 
(Figs.  1,  2a).  The  ADCP  at  M5  operated  at  600  kHz  and 
all  others  operated  at  300  kHz.  Transducer  heads  were 
approximately  0.5  m  above  the  sea  floor.  Yearlong  cur¬ 
rents  and  hydrographic  fields  were  collected  from  five 
ADCP  moorings  (M1-M5)  and  four  subsurface  string 
moorings  (S1-S4)  equipped  with  temperature  T,  con¬ 
ductivity  C,  and  pressure  P  (TCP)  sensors.  ADCPs  re¬ 
corded  nearly  full  water  column  current  profiles  every 
15  min  (an  ensemble  average  of  120  realizations)  during 
long-term  deployments.  During  the  lOP  high  temporal 
resolution,  current  and  hydrographic  data  were  collected 
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from  another  five  bottom-mounted  AD  CPs  at  stations 
M4a,  M5a,  M6,  M7,  and  M8  and  two  string  moorings  at 
S5  and  S8.  Samples  were  averaged  every  minute  at  M4a, 
M5a,  and  M6  and  every  2  min  at  M7  and  M8.  Vertical 
bin  sizes  were  1  m  at  M5a,  2  m  at  M4a,  and  4  m  at  M6, 
M7,  and  M8.  Accuracy  for  the  ADCP  measurements  is 
0.5%  of  the  observed  water  velocity  ±  0.5  cms“^. 

Barny  moorings  were  also  equipped  with  Ppods  (Mourn 
and  Nash  2008).  Ppods  are  modified  Paroscientfic  pres¬ 
sure  sensors  and  have  a  precision  of  about  0.14  mm  or 
IPa  (IPa  =  lNm“^).  A  total  of  five  Ppods  were  de¬ 
ployed  during  MORT  (Fig.  2a)  at  M2,  M4,  M5,  M7,  and 
M8,  and  their  pressure  is  referred  to  as  Pm2,  ^m4,  ^m5, 
Pm7,  and  Pms-  Three  of  the  Ppods  (Pm2,  Pm4,  and  Pms) 
were  deployed  along  a  line  over  the  bank  for  a  period  of 
nearly  6  months  (between  June  and  December  2011), 
and  the  other  two  (Pm?  and  Pms)  were  deployed  for 
2  weeks  in  June  2011  during  the  lOP.  We  noted  mis¬ 
matches  in  time  stamps  and  nonuniform  sampling  rates 
in  the  long-term  deployed  Ppods.  The  first  step  was  to 
obtain  correct  time  stamps  and  sampling  rates.  The  time 
stamps  of  Pm2  and  Pm4  were  accurate,  but  their  sampling 
rates  deviated  from  1  Hz,  indicating  drifts  in  the  sam¬ 
pling  clocks.  Note  that  the  average  sampling  rate  was 
0.97  Hz.  To  fix  sampling  rates  caused  by  clock  drifts,  we 
resampled  (i.e.,  interpolated)  Pm2  and  Pm4  data  at  1-Hz 
sampling  rates  between  accurate  time  stamps  of  starting 
and  ending  times.  We  noted  that  the  accurate  clock  in 
Pm5  was  broken  in  August  2011,  but  the  sampling  clock 
had  an  accurate  time  stamp  with  a  sampling  rate  of  1  Hz; 
when  compared  with  the  rest  of  the  Ppods,  the  sampling 
clock  in  Pm5  was  consistent  with  the  other  accurate 
clocks  in  Pm2  and  Pm4-  Therefore,  we  used  the  time 
stamp  of  the  sampling  clock  as  the  accurate  time  stamp 
for  Pm5. 

Four  TCP  string  moorings  (S1-S4)  were  deployed 
close  to  M1-M4  during  the  first  half  of  the  deployment 
and  only  three  (S2-S4)  were  redeployed  near  M2-M4 
during  the  second  half  (Fig.  2a).  Each  string  mooring 
contained  8  to  12  TCP  sensors  that  were  approximately 
equally  spaced  between  7  and  12  m  below  the  surface 
and  1  m  above  the  bottom  (Teague  et  al.  2013).  We  do 
not  address  processing  and  data  quality  issues  here, 
since  those  topics  were  discussed  in  Teague  et  al. 
(2013). 

For  the  current  data,  we  used  local  coordinates  with 
the  U  component  oriented  toward  the  axis  of  the  M2-M4 
mooring  line  and  the  V  component  oriented  perpen¬ 
dicular  to  the  mooring  line  (Fig.  2).  We  rotated  ADCP 
currents  by  19°  counterclockwise  to  the  horizontal  di¬ 
rection  to  match  with  the  direction  of  the  form  drag 
estimated  in  section  4.  Low-frequency  east-west  flow 
showed  a  two-layer  structure  around  the  bank  (Teague 


et  al.  2013).  A  similar  flow  pattern  was  found  in  U. 
Figure  3  shows  the  depth-averaged  low-frequency 
[<0.02  cycles  per  hour  (cph)]  currents  (Uu,  Vjj)  in  the 
upper  50-m  layer.  The  flow  Uu  was  nearly  eastward 
during  the  first  half  of  the  record  and  changed  direction 
after  yearday  263  (Figs.  3a,b).  Currents  in  the  lower 
layer  Vl)  were  smaller  than  in  the  upper  layer  and 
were  highly  influenced  by  the  bank  (Figs.  3c,d;  Teague 
et  al.  2013).  Mean  currents  approximately  20  m  above 
the  bottom  followed  the  bottom  topography  (Teague 
et  al.  2013).  Shipboard  and  moored  velocity  measure¬ 
ments  made  during  the  lOP  (2-13  June  2011)  showed 
generation  of  submesoscale  and  mesoscale  motions  and 
reversals  of  bottom  currents  on  the  lee  side  of  the  bank 
(Jarosz  et  al.  2014,  manuscript  submitted  to  /.  Geophys. 
Res.).  Based  on  all  the  ADCP  velocity  records,  we  noted 
that  the  currents  in  the  upper  50-m  layer  at  Ml  were  least 
affected  by  the  bank.  Ml  was  located  at  the  northwest 
corner  of  the  bank  where  the  bottom  slope  was  less  steep 
compared  to  slopes  south  and  east  side  of  the  bank 
(Fig.  2).  Near-surface  intensified  high-frequency  currents 
(>0.2  cph)  were  observed  over  the  top  of  the  bank  and 
were  closely  related  with  wind  and  surface  waves 
(Wijesekera  et  al.  2013). 

The  spectra  oi  Uu  and  Vu  show  several  distinctive 
frequency  bands  associated  with  tidal/inertial  waves  and 
subtidal  flows  (Fig.  4).  The  near-inertial  period  (25.52  h) 
lies  between  the  Ki  (23.93  h)  and  Oi  (25.  82  h)  tidal 
periods.  Teague  et  al.  (2013)  reported  that  near-inertial 
currents  were  as  large  as  20cms“^,  especially  during 
summer  months,  and  that  barotropic  diurnal  tides  were 
less  than  about  4cms“^.  The  tidal-inertial  components 
of  east-west  and  north-south  directions  had  similar 
magnitudes,  while  the  semidiurnal  components  in  the 
north-south  direction  were  stronger  than  the  east-west 
components  (Fig.  4).  The  subtidal  band  contains  4-6-, 
10-12-,  and  16-17-day  fluctuations.  There  was  no  clear 
spring-neap  tidal  signal  (period  —14  day)  in  the  subtidal 
band  (Fig.  4).  In  the  northern  GoM,  the  4-6-day  signal 
may  relate  to  synoptic-scale  weather  systems  (e.g., 
Donohue  et  al.  2006).  Donohue  et  al.  (2006,  2008)  ob¬ 
served  16-  and  14-day  oscillations  in  bottom  pressure  in 
the  north  central  GoM,  but  these  signals  were  unex¬ 
plained.  However,  Donohue  et  al.  (2008)  further  noted 
that  the  16-day  signal  is  in  phase  with  water  level  varia¬ 
tions  measured  by  coastal  tide  gauges  on  the  west  Florida 
shelf. 

Near-surface  temperature  showed  seasonal  warming 
in  summer  and  cooling  in  late  fall  (Fig.  3e).  The  advec- 
tion  of  cold  water,  especially  near  the  bottom,  appeared 
after  yearday  240  [lasting  close  to  the  end  of  the  obser¬ 
vational  period  (Fig.  3e)]  and  changed  the  vertical 
structure  of  temperature  and  density  fields,  which  in  turn 
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Fig.  3.  Time  series  of  48-h  low-pass  filtered  upper-layer  (above  50  m)  and  lower-layer  (below 
50  m)  velocities  (cms“^).  (a)  Upper-layer  velocity  along  the  direction  of  the  mooring  axis  M2- 
M4  (C/[/).  (b)  Upper-layer  velocity  perpendicular  to  the  mooring  axis  {Vu)-  (c)  Lower-layer 
velocity  along  the  direction  of  mooring  axis  M2-M4  {Ui).  (d)  Lower-layer  velocity  perpen¬ 
dicular  to  the  mooring  axis  {Vi}.  (e)  Color  image  of  depth-time  series  of  15-min  averaged 
temperature  (°C)  at  S2. 


altered  bottom  pressure  across  the  bank  on  subseasonal 
to  seasonal  time  scales.  Teague  et  al.  (2013)  suggested 
a  plausibility  of  advection  of  cold  water  by  eddies  in  the 
region.  In  the  following  analysis,  we  focus  on  motions 
that  can  be  well  resolved  by  our  data,  and  therefore  we 
limit  the  computation  of  the  form  drag  to  tidal  (12-25  h) 
and  subtidal  (4-17  days)  motions. 

b.  Ship  surveys:  Currents,  hydrography,  and 
microstructure 

Ship  surveys  were  conducted  from  the  National  Oce¬ 
anic  and  Atmospheric  Administration  (NO  A  A)  R/V 


Manta,  a  25-m  catamaran  vessel  out  of  Galveston,  Texas. 
A  300-kHz  RDI  AD  CP  was  mounted  on  the  side  of  the 
boat  and  recorded  currents  with  a  vertical  bin  size  of  4  m. 
Velocities  were  processed  using  the  University  of  Hawaii’s 
Common  Oceanographic  Data  Access  System  (CODAS; 
Firing  et  al.  1995)  to  correct  for  ship  motion  and  were 
averaged  over  4-min  intervals.  A  BioSonics  DT-X  digital 
scientific  echosounder  was  used  to  acquire  bathymetric 
data  based  on  acoustic  backscatter.  The  echosounder 
sampled  at  a  rate  of  3  Hz  with  a  frequency  of  123  Hz. 
The  average  horizontal  sampling  distance  was  approxi¬ 
mately  25  cm. 
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Fig.  4.  Depth-averaged  frequency-weighted  velocity  spectra  in 
the  upper  50  m  based  on  time  series  of  velocity  shown  in  Figs.  3a 
and  3b.  (a)  t/[/  and  (b)  Vu-  Spectra  at  Ml,  M4,  and  M5  are  plotted. 
The  vertical  lines  denote  frequencies  of  M2,  Ki,  Oi,  /  (inertial), 
4-day,  and  17-day  motions. 


Additionally,  vertical  profiles  of  small-scale  velocity 
and  temperature  along  with  finescale  T,  C,  and  P  fields 
were  collected  from  a  free-falling  Rockland  Scientific 
vertical  microstructure  profiler  (VMP),  equipped  with  a 
pressure  sensor,  shear  probes,  a  fast  thermistor,  micro¬ 
conductivity,  and  Sea-Bird  Electronics  temperature  and 
conductivity  sensors  (e.g.,  Wolk  et  al.  2002).  These 
profiles  were  collected  while  the  R/V  Manta  was  drifting 
or  steaming  slowly  along  quasi-straight  lines  over  the 
bank  at  speeds  of  0.5  to  1  m  s“^  (blue  and  green  lines  in 
Fig.  2a).  The  velocity  shear  accuracy  was  5%,  and  the 
VMP  probe  had  an  average  descent  speed  of  0.8 ms“^, 
producing  shear  measurements  approximately  every 
0.2  cm.  Shear  measurements  were  utilized  to  estimate 
turbulent  kinetic  energy  (TKE)  dissipation  rates  using 
8  =  1 .5v{dul dzf ,  where  v  is  the  kinematic  viscosity  of 
water  estimated  from  concurrent  temperature  observa¬ 
tions,  and  duldz  is  the  velocity  shear  (Gregg  1987).  The 
velocity  shear  variance  was  calculated  by  integrating  the 
spectrum  for  the  wavenumber  range  from  2  to  30  cycles 
per  meter  (cpm),  while  the  spectrum  was  estimated  in 
bins  of  1024  shear  data  records,  with  a  512  record 
overlap.  We  noted  that  more  than  90%  of  the  spectral 
variance  was  between  2  and  30  cpm  wavenumbers.  In 
this  study,  we  focus  on  the  dissipation  estimates  near 
the  bottom.  A  detailed  discussion  of  dissipation  esti¬ 
mates  and  their  distribution  during  the  lOP  is  given  in 
Jarosz  et  al.  (2014,  manuscript  submitted  to/.  Geophys. 
Res.). 


3.  Form  and  frictional  drags 

a.  Form  drag 

The  form  drag  is  defined  as  the  spatial  integral  of  the 
product  of  the  bottom  pressure  anomaly  and  the  lateral 
slope  of  the  bathymetry  (e.g..  Gill  1982,  p.  145;  Baines 
1995).  The  form  drag  along  a  line  of  pressure  measure¬ 
ments  is  (e.g.,  McCabe  et  al.  2006;  Warner  et  al.  2013) 


^form(0- 


'^2 

p'j^{x,y,t)h^dx  [Nm“^],  (1) 

Jx^ 


where  p'^  is  the  bottom  pressure  anomaly  in  N  m~^,  /z  is 
the  bottom  height  in  meters,  is  the  lateral  gradient 
of  the  bank  in  the  direction  where  the  drag  is  estimated,  t 
is  the  time,  x  and  y  are  the  orthogonal  coordinates  over 
the  bank,  and  Xi  and  X2  are  points  of  the  equal  depth  on 
either  side  of  bank  (Fig.  2).  The  power  (Pform)  associated 
with  Oform,  or  the  rate  of  energy  loss  of  the  background 
flow  due  to  the  form  drag,  is  defined  as  the  product  of 
/)form  and  the  free-stream  velocity  undisturbed  by  the 
topographic  feature  Um  (Gill  1982,  section  8.7;  Warner 
et  al.  2013),  where 


P,  =  Dr  U  . 

form  form  m 


(2) 


The  drag  coefficient  Cdf  associated  with  the  form  drag 
[(!)]  is  defined  as  (e.g.,  Lott  and  Miller  1997;  Vosper 
et  al.  1999;  Warner  and  MacCready  2009) 


r  =  — 


^form^^O 


O.SpJU  \U  ’ 

^01  ml  m 


(3) 


where  pq  is  the  seawater  density,  and  ho  is  the  topo¬ 
graphic  height.  Here  a  minus  sign  is  included  since  the 
drag  and  the  velocity  have  opposite  signs.  The  term  Cdf 
is  an  0(1)  quantity  and  is  referred  to  as  the  bluff-body 
drag  coefficient  (Warner  et  al.  2013).  The  quantity  Cdf 
is  widely  used  in  engineering  applications.  To  give  an 
oceanographic  context,  we  also  defined  a  bulk  (form) 
drag  coefficient  associated  with  Oform  as  an  equiv¬ 
alent  frictional  drag  coefficient  based  on  mean  flow  ve¬ 
locity  Um  (e.g.,  Warner  et  al.  2013),  where 


^form^^Q 


-O.SpJU  \u 

m\  rr, 


(4) 


and  lo  is  the  topographic  length.  Equation  (4)  defines  the 
bulk  drag  coefficient  in  the  same  way  as  the  frictional 
drag  coefficient,  although  the  near-bottom  velocity  is 
typically  used  in  computing  the  frictional  drag  resulting 
from  bottom  roughness  elements  [see  (13)].  The  form 
drag  coefficient  Cd  can  be  treated  as  an  enhanced 
frictional  drag  coefficient  from  a  numerical  modeling 
perspective.  From  (3)  and  (4),  Cdf/^e)  =  hlho,  and  for 
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the  EFGB  (Fig.  2)  with  ~  50  m,  and  /q  ~  9  km,  Cop/ 
Cd  ~  180. 

The  first  step  in  estimating  the  form  drag  is  to 
determine  p'^  from  bottom  pressure  measurements. 
The  bottom  pressure  sensors,  deployed  during  MORT, 
measured  the  total  pressure  at  the  seafloor  pb,  which 
includes  atmospheric  pressure,  pressure  resulting  from 
surface  elevation  and  water  density,  and  non¬ 
hydrostatic  pressure.  Following  Mourn  and  Smyth 
(2006)  and  Warner  et  al.  (2013),  we  express  the  total 
pressure  p(X,  t)  at  depth  z  by  integrating  the  vertical 
momentum  balance  from  z  to  the  sea  surface  height 
T?(x,  y,  t): 


p{X,  t)  =  +  gPQ[vix,y,  0  -  2]  +  ^  J  pW,  0  dz' 

•17  p 

+  g  p(X,  t)  dz'  +  Pq  DwIDt  dz' ,  (5) 

Jz  Jz 


where  p(X,  t)  varies  in  all  directions;  X  =  (x,  y,  z)  rep¬ 
resents  a  right-handed  coordinate  system  with  the  x  axis 
parallel  to  the  bottom  pressure  sensor  array,  the  y  axis 
perpendicular  to  the  sensor  array,  the  z  axis  positive 


upward,  and  z  =  0  at  the  bottom;  pAtm  is  the  atmospheric 
pressure  at  the  surface  and  is  assumed  to  be  uniform 
over  the  domain;  D/Dt  is  the  material  derivative;  w  is  the 
vertical  velocity;  and  g  is  the  gravitational  acceleration 
(9.81  ms“^).  The  density  is  expressed  as  p(X,t)  = 
Pq  +  p(z,  t)  +  p'(X,  t),  where  po  =  constant,  p(z,  t)  is  the 
background  stratification,  and  p'(X,  t)  is  the  density 
fluctuation.  The  sea  surface  can  be  divided  into  four 
components: 

'n{x,y,  t)  =  rj(t)  +  VTf(^,y,  t)  +  VTM,y,  0  +  'n'{x,y,  t) , 

(6) 

where  p(t)  is  the  spatially  averaged  sea  surface  height; 

0  is  the  contribution  to  sea  surface  height 
associated  with  the  sea  surface  slope  resulting  from  rota¬ 
tional  effects.  Here  barotropic  pressure  gradients  associ¬ 
ated  with  surface  elevation  are  balanced  by  the  Coriolis 
force;  17^7 (x,  y,  t)  is  the  surface  elevation  resulting  from 
accelerations/decelerations  of  flow  over  the  bank,  and 
T7'(x,  y,  t)  is  the  remaining  sea  surface  height  anomaly 
(e.g.,  Warner  et  al.  2013).  By  combining  (5)  and  (6)  we 
obtain 


p{X,  t)  =  +  gPoT?(f)  +  gPoVj-fix, y,  t)  +  gpQVj.i(x, y,  t)  +  gp^vix, y,  t)  -  gp^z 

+  ^J  p{z',t)dz  +  g^  p{X,t)dz'  +  PQ^  Dw/Dtdz'.  (7) 


Pressure  terms,  which  vary  in  the  x  direction,  contribute 
to  the  form  drag  given  in  (1).  The  pressure  anomaly  that  is 
dynamically  relevant  to  the  flow  over  the  bank  is  given  by 

p'(X,  t)  =  gPQ'r)ji(x,y,  t)  +  gpQT]'(x,y,  t) 

+  g  p'iX,  t)  dz'  +  Po  ?  DwIDt  dz' .  (8) 

Jz  Jz 

Rotational  contributions  to  sea  surface  height  are  ex¬ 
cluded  because  their  associated  velocity  component  is 
orthogonal  to  the  direction  of  the  form  drag. 

b.  Frictional  drag 

Frictional  stress  resulting  from  turbulence  and  viscous 
effects  (Schlichting  1962;  Tennekes  and  Lumley  1972)  is 

Pm=Pi^^)~P^^  [Nm"2],  (9) 

where  U  is  the  mean  lateral  velocity,  u  and  w  are  the 
lateral  and  vertical  velocity  fluctuations,  v  is  the  kine¬ 
matic  viscosity,  and  z  is  the  distance  from  the  sea  floor 
(positive  upward)  as  defined  above.  The  angle  brackets 


denote  the  average  over  turbulent  motions.  In  the 
presence  of  strong  surface  waves,  the  turbulent  bottom 
stress  —p{uw)  over  the  continental  shelf  and  shallow 
banks  may  include  a  stress  component  from  surface 
wave-induced  oscillatory  currents  (e.g..  Grant  and 
Madsen  1982).  The  surface  wave-driven  stress  is  sig¬ 
nificantly  larger  over  the  shelf  during  severe  storms  (e.g., 
Wijesekera  et  al.  2010),  but  in  the  following  analysis,  we 
limit  our  bottom  drag  computations  to  mean  currents, 
since  hurricane  strength  winds  were  not  present  over  the 
FFGB  during  the  observational  period.  As  discussed 
below,  we  estimated  bottom  roughness  lengths  and 
frictional  drag  coefficients  when  microstructure  obser¬ 
vations  were  available  for  10  days  in  June  and  August 
2011.  Winds  were  less  than  6ms“^  and  surface  wave 
heights  were  0.2-0.5  m  (Wijesekera  et  al.  2013).  Surface 
wave-induced  bottom  currents  were  negligible  over  the 
top  of  the  bank.  Therefore,  the  surface  wave-induced 
bottom  frictional  effects  or  “the  apparent  roughness 
lengths”  were  not  a  significant  factor.  The  wall  boundary 
layer  theory  provides  an  estimate  of  the  roughness 
length  of  bottom  roughness  elements. 
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Sufficiently  far  from  the  boundary  z  >  vlu  or  at 
sufficiently  large  Reynolds  numbers  uzlv  >  1,  a  layer 
(referred  as  the  inertial  layer)  exists  in  which  the 
overall  dynamics  of  turbulence  is  independent  of  vis¬ 
cosity.  In  this  layer,  =  p{uw)  =  —  r,  where  r  is  the 
bottom  stress.  The  fiux  term  can  be  rewritten  as 
a  downgradient  fiux,  Fm  =  -pK^dUldz,  where  Km  is 
the  turbulent  transfer  coefficient  having  dimensions  of 
velocity  and  length  and  dU/dz  is  the  vertical  gradient  of 
the  mean  fiow.  For  the  inertial  layer.  Km  =  ku^z, 
|t|  =  and 


where  k  is  von  Karman’s  constant  (0.41),  and  is  the 
water  friction  velocity  that  is  independent  of  z.  The 
turbulent  fluxes  in  the  inertial  layer  are  nearly  inde¬ 
pendent  (within  10%-20%)  of  the  distance  from  the 
boundary.  By  integrating  (10),  we  can  obtain  a  mean 
velocity  profile: 

[/  =  -log(E),  (11) 

K  V^O/ 

where  Zo  is  the  roughness  length  and  is  determined  by 
the  geometry  of  the  surface.  The  term  zo  is  a  small  frac¬ 
tion  of  the  frictional  element  heights;  for  example,  for 
closely  packed  sand  grains  of  size  A,  zo  =  A/30  (Nikuradse 
1950). 

Suppose  we  have  observations  in  the  inertial  layer  at 
z  =  Zfl,  where  U(za)  =  Ua,  then  squaring  (11),  we  solve 
for  the  magnitude  of  the  bottom  stress: 

\T\=p[KHog(zJZo)~^]Ul  =  pC,Ul,  (12) 

where 

Ca  =  [Kl\og{zJZ(^)f  (13) 

is  the  frictional  drag  coefficient.  Within  the  inertial  layer 
the  production  of  TKE  by  sheared  currents  approxi¬ 
mately  balances  the  dissipation  rate,  and  thus 

=  {sKz)^'^ ,  (14) 

where  e  is  the  TKE  dissipation  rate  (e.g.,  Dewey  and 
Crawford  1988).  The  roughness  length  becomes 

Zq  =  (15) 

In  the  following,  we  use  velocity  and  TKE  dissipation 

rate  profiles  to  compute  the  roughness  length  and  the 
associated  frictional  drag  coefficient  over  the  EFGB. 


Fig.  5.  (a)  Demeaned  bottom  pressure  {P  -  {P))  at  M2,  M4,  M5, 
M7,  and  M8  between  2  Jun  (yearday  152)  and  13  Jun  (yearday  163) 
2011,  where  (P)  is  the  time-averaged  pressure.  Pressure  at  M7 
(south  of  the  bank;  Fig.  2a),  marked  in  cyan  differs  from  the  rest  of 
the  observations,  (b)  Depth-averaged  velocity  in  the  direction  of 
the  mooring  axis  {Ud)  and  (c)  velocity  perpendicular  to  the 
mooring  axis  {Vd)  at  M1-M5. 

4.  Evaluation  of  form  drag 

a.  Bottom  pressure  anomaly  from  Ppods 

Figure  5  shows  bottom  pressure  fluctuations  at  five 
mooring  sites  (Figs.  1,  2)  between  2  June  (yearday  152) 
and  13  June  (yearday  163)  2011.  The  pressure  fluctua¬ 
tions  were  constructed  by  removing  the  temporal  mean 
for  a  given  record  length  of  12  days.  The  demeaned 
pressure  fluctuations  show  diurnal  tides  and  a  segment 
of  the  spring-neap  cycle.  Tides  around  and  over  the  bank 
had  similar  phases  (Fig.  5a)  since  the  spatial  scales  of  tidal 
motions  are  much  larger  than  the  lateral  spacings  be¬ 
tween  moorings.  Unlike  pressure,  the  depth- averaged 
velocities  at  M1-M5  were  not  always  in  phase,  especially 
at  the  beginning  of  the  records  (Figs.  5b,c).  While  pres¬ 
sure  fluctuations  were  dominated  by  diurnal  tides,  ve¬ 
locity  fluctuations  contained  combinations  of  tides  and 
near-inertial  waves. 

Pressure  spectra  calculated  at  M2,  M4,  and  M5  (Fig.  6) 
show  several  distinctive  features  including  the  surface- 
gravity  wave  band  between  0.08  and  0.2  Hz,  the  infragravity 
wave  band  between  0.002  and  0.08  Hz,  the  semidiurnal 
M2  tide  (period  =  12.42  h),  and  the  diurnal  tidal  bands  Ki 
(period  =  23.93  h)  and  Oi  (period  =  25.82  h).  Our  pres¬ 
sure  spectra  also  indicate  variability  in  the  synoptic  time 
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Fig.  6.  (a)  Frequency  spectra  of  bottom  pressure  at  M2,  M4,  and  M5.  The  spectra  were  based 
on  the  entire  pressure  record  with  a  sampling  rate  of  1  Hz.  Spectra  were  averaged  into 
frequency  bins,  (b)  Enlarged  plot  showing  diurnal  tides  Oi  and  Ki  and  semidiurnal  tide  M2. 


scale  (4-6  days)  and  in  10-17-day  fluctuations,  similar  to 
velocity  spectra  shown  in  Fig.  4.  The  amplitude  of  a  given 
tidal  frequency  was  estimated  by  integrating  the  pressure 
spectrum  for  a  specified  spectral  band,  and  the  ampli¬ 
tudes  of  Oi,  Ki,  and  M2  tides  are  0.17,  0.17,  and  0.06  m, 
respectively.  The  pressure  spectrum  (Fig.  6)  has  a  slope  of 
about  -1.5  between  10“"^  and  10“^  Hz.  Teague  et  al. 
(2013)  examined  velocity  components  of  tides  and  near- 
inertial  waves  (period  =  25.52  h).  They  reported  that  tidal 
velocities  were  smaller  than  near-inertial  velocities,  and 
magnitudes  of  Ki  and  Oi  were  about  3-5cms“^.  The 
magnitude  of  the  M2  tidal  velocity  was  about  2.5  cms~^. 

We  adopted  the  following  procedure  to  compute 
pressure  anomalies.  First  high-frequency  pressure  fluc¬ 
tuations  and  spikes  were  removed  by  filtering  raw  pres¬ 
sure  to  remove  frequencies  higher  than  10“^  Hz.  Second, 
seasonal  and  subseasonal  motions  with  time  scales 
greater  than  30  days  were  filtered  out,  since  these  motions 
are  not  well  resolved,  and  our  focus  is  on  tidal  and  sub- 
tidal  (4-17  day)  oscillations.  These  bandpass-filtered 
pressure  fluctuations  at  M2,  M4,  and  M5  are  referred  as 
Pmi,  Pma,  and  Pm5*  The  30-day  high-pass  filter  basically 
removed  the  largest  term  (pogz)  in  (7).  The  sea  surface 
elevation  PQg^(t)  was  removed  by  subtracting  the  refer¬ 
ence  or  the  background  pressure  at  a  location  away  from 
the  bank.  Here  we  used  Pmi  as  the  reference  pressure 
record.  This  choice  of  the  reference  pressure  does  not 
affect  the  form  drag  calculation,  and  any  Ppod  pressure 


could  have  been  used  as  a  background.  The  value  Pog^(t) 
is  spatially  uniform  and  therefore  does  not  affect  the  form 
drag  calculation.  Pressure  anomalies  at  M4  (APm4)  and  at 
M5  (APms)  were  obtained  by  removing  Pm2  from  Pm4 
and  Pm5,  where  APm4  =  Pm4  -  Pm2  and  APm5  =  ^m5  “ 
Pm2  (Fig*  7).  Note  the  pressure  anomaly  at  M2,  where 
APm2  became  zero.  The  subtraction  of  the  reference 
pressure  also  removed  the  atmospheric  pressure  pAtm, 
since  pAtm  was  assumed  to  be  uniform  over  the  bank.  The 
background  stratification  term  g  j^p(P,  t)dz'  was  not 
computed  because  the  time-dependent,  depth- averaged 
density  was  not  well  known.  Time-depth  series  of  hy¬ 
drographic  data  from  three  TCP  moorings  at  S2,  S3,  and 
S4  were  available,  but  lateral  and  vertical  resolutions  of 
the  density  field  were  very  coarse  to  make  an  accurate 
estimate  of  the  background  stratification  term.  This  term 
is  spatially  uniform,  and  therefore  the  subtraction  of  the 
reference  pressure  Pm2  removes  the  term  g  p(z',  t)  dz' 
from  pressure  anomalies  at  M4  and  M5.  The  sea  surface 
tilt  due  to  rotational  effects  was  estimated  from  the 
geostrophic  balance  gd^^j-ldx  =  fV,  where  V  is  the  low- 
pass  filtered  (periods  greater  than  48  h),  depth- averaged 
velocity  perpendicular  to  the  mooring  line,  and  /  is  the 
Coriolis  parameter.  We  noted  that  V  was  nearly  uniform 
around  the  bank  (not  shown).  Therefore,  the  pressure  in¬ 
duced  by  the  sea  surface  height  due  to  rotational  effects  at 
a  given  location  along  the  mooring  line  was  approximated 
from  gPoTJTj  ^  po/C(x  —  xm2)*  For  /  =  6.84  X  10”^  s~^. 
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Fig.  7.  Bottom  pressure  anomalies  (thin  gray  lines)  at  (a)  M4, 
APm4  =  Pma  -  Pmi,  and  (b)  M5,  APms  =  ^m5  -  Pmi-  Barotropic 
pressure  anomaly  estimated  from  rotational  effects,  gpo'^lr/ 
Pof^i^  ~  tm2),  are  marked  by  thick  black  lines. 

y  =  0.1  m  pq  =  1000  kg  m“^,  and  x  =  Xm4  with  (xm4  ~ 
^Mi)  =  9  km,  the  resulting  pressure  anomaly  at  M4  is  as 
large  as  100  Nm“^  (Fig.  7a). 

As  described  above,  the  bottom  pressure  anomalies 
west,  east,  and  in  the  middle  of  the  bank  (M2,  M4,  and 
M5,  respectively)  were  obtained  by  (i)  bandpass  filtering 
of  pressure  between  1000  s  and  30  days,  (ii)  subtracting  a 
selected  reference  pressure  outside  the  bank  (i.e.,  pres¬ 
sure  at  M2),  and  (iii)  subtracting  the  pressure  due  to 
rotational  effects.  Bottom  pressure  anomalies  at  M2, 
M4^and  M5  are,  respectively,  ^  0,  =  APm4- 

PofV(xM4  -  xmi),  andp^j^3  =  APm5  -  PofV  (xms  - 
The  bottom  pressure  anomalies  were  dominated  by  the 
combined  external,  internal,  and  inertial  pressure 
anomalies  (APm4  and  APms)-  The  rotational  effect  that 
depends  on  the  relative  distance  from  M2  is  important, 
especially  at  M4,  but  it  is  a  smaller  term  at  M5  (Fig.  7). 
Bottom  pressure  fluctuations  (Fig.  7)  show  multiple  time- 
scale  variability  in  APm4  and  APm5-  They  were  about  50- 
100  Nm“^  at  tidal  frequencies  and  about  50-200  Nm“^ 
for  low-frequency  motions  (periods  ~  4-17  days). 

b.  Form  drag 

We  used  three  pressure  sensors  at  M2,  M5,  and  M4 
(Fig.  2a)  to  compute  the  form  drag  over  the  bank.  The  drag 
was  computed  along  the  mooring  line  that  was  oriented 
about  19°  counterclockwise  from  the  east-west  direction 
(Fig.  2).  The  topographic  section  along  the  mooring  line  is 
close  to  “top  hat”-shaped  bathymetry.  In  general,  the 


EFGB  has  steeper  slopes  along  the  east  and  south 
boundaries  and  a  significant  portion  of  the  bank  top  is 
relatively  flat,  except  near  the  peak  (Fig.  2).  The  major 
challenge  is  to  evaluate  the  form  drag  [(1)]  from  only 
three  pressure  measurements.  Hafner  and  Smith  (1985) 
estimated  the  form  drag  over  the  European  Alps  using  a 
limited  number  of  pressure  measurements.  They  trans¬ 
formed  the  surface  integral  [(1)]  to  the  flux  integral  and 
evaluated  the  form  drag  as  the  product  of  the  horizontal, 
locally  constant  pressure  gradient  and  the  volume  of  the 
mountain  terrain.  This  is  the  equivalent  of  assuming  that 
the  local  pressure  gradient  is  linear  between  two  nearby 
observation  sites.  Hafner  and  Smith  (1985)  discussed  the 
uncertainty  of  pressure  that  was  determined  by  a  linear 
formula  based  on  their  surface  stations. 

The  form  drag  is  large  in  regions  where  both  bathy¬ 
metric  slopes  and  pressure  anomalies  are  significant. 
Warner  et  al.  (2013)  examined  the  sensitivity  of  a  num¬ 
ber  of  pressure  measurements  by  evaluating  the  integral 
[(1)]  over  a  sloping  headland,  TTP  in  Puget  Sound, 
Washington,  and  found  that  four  sensors  (two  on  either 
slope  of  the  ridge;  Fig.  Al)  can  provide  an  accurate  es¬ 
timate.  We  also  examined  sensitivity  for  our  three-point 
measurement  setup  (one  on  the  top  and  one  on  either 
side  of  the  ridge;  Fig.  Al)  using  Warner’s  (2012)  numer¬ 
ically  generated  pressure  fields  over  the  TTP  ridge.  As 
described  by  Warner  (2012),  at  TTP,  the  large  bottom 
pressure  anomalies  were  found  where  the  topography  was 
the  steep.  The  model  was  able  to  capture  some,  but  not  all 
of  these  pressure  anomalies.  The  spatial  resolution  of  the 
model  and  the  required  smoothing  of  the  topography  for 
the  model  to  run  were  major  factors  for  simulating  an 
accurate  pressure  field  over  the  bank  (Warner  2012). 
Since  the  Ppods  in  MORT  were  not  located  on  the  sloping 
sections  of  the  topography,  some  of  the  largest  pressure 
anomalies  could  have  been  missed. 

A  sensitivity  analysis  described  in  the  appendix  allows 
the  examination  of  uncertainties  in  our  computation 
based  on  three  pressure  measurements,  although  the 
geometrical  and  dynamical  settings  of  this  study  and 
Warner  et  al.’s  (2013)  are  not  identical  and  the  model 
results  may  not  perfectly  capture  the  distribution  of 
bottom  pressure  anomalies  due  to  numerical  limitations. 
We  compared  form  drag  estimates  based  on  multiple 
pressure  measurements  versus  three  pressure  measure¬ 
ments  as  described  in  the  appendix.  We  used  three 
methods  to  compute  the  form  drag  over  the  bank. 
Method  I  computes  the  form  drag  by  approximating 
bottom  pressure  for  a  given  slope  as  the  mean  pressure 
between  the  top  and  the  outer  edge  of  the  ridge  and  then 
multiplying  it  by  the  mean  slope.  Method  II  computes 
the  form  drag  as  a  product  of  the  observed  slope  and  a 
linearly  interpolated  pressure  between  the  outer  edges 
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Fig.  8.  Form  drag  estimated  from  the  three  methods  described  in 
the  appendix, 

and  the  top  of  the  ridge.  Method  III  computes  the  form 
drag  as  a  product  of  the  mean  pressure  for  a  given  slope 
and  the  maximum  magnitude  of  that  slope.  The  estimates 
from  methods  I  and  II  were  comparable,  but  their  max¬ 
imum  magnitudes  were  a  factor  of  2  smaller  than  those 
obtained  from  the  method  III  and  the  fully  resolved  so¬ 
lution.  Here  we  applied  all  three  methods  to  compute  the 
form  drag  along  the  M2-M4  mooring  line,  and  results  are 
illustrated  in  Fig.  8.  Topographic  gradients  of  the  bank 
are  small  or  near  zero  over  the  top  and  outside  the  100-m 
bathymetry  contour  compared  to  the  two  major  slopes  of 
the  bank  (Figs,  la,  2).  Therefore,  we  can  expect  major 
contributions  to  the  form  drag  from  sloping  boundaries. 
The  EFGB  is  asymmetric,  especially  toward  the  southern 
end,  where  the  lateral  slopes  over  the  bank  are  not  neg¬ 
ligible.  We  may  not  compute  the  form  drag  completely 
since  our  computation  is  limited  to  the  east-west  di¬ 
rection.  Drags  estimated  from  linearly  interpolated  pres¬ 
sure  and  from  mean  pressure  are  comparable.  The  form 
drag  estimated  from  method  III  is  twice  large  as  those  that 
calculated  from  methods  I  and  II,  similar  to  the  findings  in 
the  appendix,  indicating  the  sensitivity  to  the  slope. 

In  the  following  analyses,  we  used  the  form  drag 
computed  from  method  II,  similar  to  Hafner  and  Smith 
(1985),  even  though  we  may  underestimate  the  expected 
value  by  a  factor  of  2.  To  examine  the  time  variability  of 
the  total  form  drag,  we  examine  tidal  and  subtidal 
components  (time  scales  greater  than  2  days)  separately. 
The  corresponding  drag  estimates  based  on  method  II 
are  shown  in  Figs.  9a  and  9b.  The  tidal  band  contains 
diurnal  and  semidiurnal  tides,  and  the  subtidal  band 
contains  multiple  time  scales  (4-17  days).  Here  it  is 


assumed  that  form  drag  can  be  expressed  as  a  linear 
summation  of  several  components.  However,  this  as¬ 
sumption  may  have  limitations  when  the  processes  that 
generate  form  drag  become  highly  nonlinear.  Frequency 
spectra  of  Dform  were  also  computed  to  examine  multi¬ 
ple  time-scale  processes  related  to  the  total  form  drag 
(Fig.  10).  The  semidiurnal  and  diurnal  tidal  components, 
4-6-day  synoptic  time-scale  variability,  and  10-17-day 
low-frequency  variability  are  major  contributors  to  the 
total  form  drag,  similar  to  the  velocity  spectrum  shown 
in  Fig.  3.  The  most  dominant  components  of  Dform  are 
the  diurnal  and  4-17-day  band  with  drag  magnitudes  in 
excess  of  2000Nm“^  (Table  1;  Figs.  9).  Note  that  the 
mean  form  drag  for  the  tidal  band  was  computed  by 
averaging  maximum  tidal  magnitudes,  and  the  mean  for 
the  4-17-day  band  was  computed  by  averaging  maxi¬ 
mum  magnitudes  in  every  4-day  time  interval. 

The  total  form  drag  [(1),  (8)]  may  contain  inertial 
pressure  exerted  by  the  oscillatory  nature  of  the  flow 
(e.g.,  Warner  and  MacCready  2009;  Warner  et  al.  2013). 
Because  of  tidal  acceleration,  such  as  during  slack  tides, 
the  barotropic  pressure  gradient  (i.e.,  sea  surface  slope 
dr]jildx)  is  in  the  direction  of  the  tidal  acceleration 
{dUldt),  which  in  turn  generates  an  “inertial  drag”  or  an 
“apparent  form  drag”  without  flow  separation.  We  noted 
that  the  tendency  term  (dUldt)  of  the  depth- averaged 
velocity  was  nearly  uniform  over  and  around  the  bank 
(not  shown).  Therefore,  the  pressure  due  to  acceleration 
of  the  oscillatory  flow  is  approximated  as  gp^'inji  ^ 
PQdUldt(x  —  xmi),  where  U  is  the  depth- averaged  veloc¬ 
ity  along  the  M2-M4  line,  and  the  corresponding  inertial 
drags  (Diner)  for  tidal  and  subtidal  bands  are  plotted 
in  Fig.  9.  For  tidal  motions.  Diner  is  a  significant  part  of 
the  total  form  drag  (Fig.  9a),  but  for  subtidal  motions 
(Fig.  9b),  Diner  is  negligible. 

c.  Power: 

Both  currents  and  the  estimated  form  drag  consist  of 
nonstationary  processes,  and  conventional  Fourier 
cospectrum  analysis  does  not  provide  information  about 
power  variations  within  data  records.  Therefore,  the 
power  per  unit  length  Pform  is  estimated  as  a  product  of 
Dform  and  11^  [(2)],  where  Um  is  the  free-stream  velocity 
in  the  direction  of  the  form  drag  at  an  upstream  location 
where  the  flow  is  undisturbed  by  the  topographic  fea¬ 
ture.  Unfortunately,  our  observations  outside  the  bank 
were  close  to  the  walls  of  the  bank.  The  flow  below  the 
height  of  the  bank  was  impacted  by  the  topography,  and 
among  those  observations,  the  flow  in  the  upper  50  m  at 
Ml  was  the  least  affected.  Therefore,  we  used  the  upper 
50-m  averaged  velocity  D^/  at  Ml  as  the  undisturbed 
mean  flow  D^.  A  time  series  of  Pform  foi"  the  subtidal 
band  at  M1-M5  is  shown  in  Fig.  11.  For  most  occasions. 
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Fig.  9.  Total  form  drag  Dform  based  on  interpolated  pressure 
fluctuations  (method  II).  (a)  Tidal  band  and  (b)  sub  tidal  bands 
(period  >  2  days).  The  blue  line  in  (a)  and  (b)  is  the  inertial 
component  of  the  form  drag  Diner  due  to  acceleration/deceleration 
of  the  flow. 


10^  10*^  10*^  10'^  10° 


Frequency,  to  (cph) 

Fig.  10.  Variance-preserving  power  spectra  of  the  form  drag 
[(N  m“^)^  X  10~^].  Thin  lines  with  crosses  denote  95%  confidence 
limits.  Vertical  lines  denote  periods  of  Ki,  Oi,  4-day,  and  17-day 
motions. 


—  Uu  and  Zlform  were  in  phase,  indicating  kinetic  energy 
loss  in  mean  currents  due  to  the  form  drag  over  the  bank 
with  negative  Pform*  Positive  values  of  Pform  laiay  reflect 
phase  differences  between  the  drag  and  the  velocity,  but 
some  differences  may  be  attributed  to  uncertainties  in 
P^form-  The  energy  loss  (Pform  <  0)  has  a  wide  range  of 
variability  with  values  as  large  as  800  Wm“^  (Fig.  11). 
The  6-month  averaged  energy  loss  (Pform  <  0)  at  Ml  was 
about  85  Wm“^  with  95%  confidence  limits  between  76 
and  93  Wm“^  (Table  1).  If  the  energy  loss  is  computed 
from  the  depth-averaged  velocity,  the  resulting  6-month 
averaged  energy  loss  at  Ml  is  73  W  m“^,  about  14%  less 
than  the  estimate  based  on  Uu  (Table  1).  The  power 
associated  with  the  diurnal  tidal  components  was  not 
computed,  since  tidal  velocities  were  heavily  contami¬ 
nated  by  inertial  motions.  The  magnitude  of  average 
energy  loss  per  unit  area,  Pform/^o,  at  subtidal  (4  to  17 
days)  bands  is  about  10~^Wm“^.  Energy  loss  over  the 
bank  can  be  expressed  as  (Pform/^o)2lo5  where  is  the 
projection  of  the  bank  to  the  horizontal  plane.  The  long¬ 
term  average  energy  losses  for  the  4-17-day  band  is 
about  0.32  MW,  where  /q  =  9  km,  andylo  (=40.76  km^) 
is  the  projected  bank  area  encompassed  by  the  100-m 
bathymetric  contour. 

d.  Drag  coefficient 

The  bluff-body  drag  coefficient  Cqf  [(3)]  and  the  bulk 
drag  coefficient  Cu>  [(4)]  were  computed,  and  the  esti¬ 
mates  of  Cu>  are  presented  here  since  Cu>yICd  =  V^o- 


The  bulk  drag  coefficient  was  calculated  from  the  ob¬ 
served  relation  between  P/)(=— 0.5po|f/m|f>^m)  and 
P)form/^o  for  tho  sub tidal  band,  where  Um  was  set  to  f/j/  at 
Ml  (Fig.  3a).  The  relationships  between  P>form/^o  and  Fu 
were  explored  for  subtidal  flows  at  all  five  mooring  sites, 
and  the  bin-averaged  P>form/^o  is  plotted  as  a  function  of 
Fu>  for  Uu  values  at  M1-M5  (Fig.  12).  As  discussed 
earlier,  Uu  Ml  is  the  best  approximation  for  an  un¬ 
disturbed  velocity  for  estimating  C^),  but  observations  at 
M2-M5  also  produced  similar  results.  The  slope  of  the 
scatterplots  represents  the  bulk  drag  coefficient.  The 
average  Cu  based  on  a  least  squares  fit  for  If)  >  Fu> 
-10  is  about  26  X  10  .  Individual  estimates  of  Cu  (i.e., 
the  ratio  P>form/^o  Fd)  vary  between  5  X  10”^  and  180  X 
10“^  with  a  mean,  standard  deviation,  and  standard  er¬ 
ror  of  about  36  X  10“^,  20  X  10“^,  and  5  X  10“^,  re¬ 
spectively.  Here  we  excluded  a  few  data  points  for  F^  < 
—  10  at  Ml.  For  the  tidal  band  we  used  (4)  to  evaluate  Cu 
from  the  tidal  drag  (Fig.  9;  Table  1)  and  mean  tidal 
currents  from  Teague  et  al.  (2013;  Table  1);  Cu  was  not 
computed  as  a  function  of  time  because  it  is  difficult  to 
generate  a  time  series  of  accurate  tidal  flow  without 
contamination  by  inertial  currents.  The  estimated  bulk 
drag  coefficients  for  the  tidal  and  low-frequency  bands 
are  about  (26-180)  X  10“^  (Table  1). 

e.  Wave  drag  from  linear  theory 

The  approximation  P>form/fe  ^  —  0.5poCz)|l/c/|l/c/  (Fig. 
12),  especially  at  negative  values  of  Uu,  suggests  that  the 
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Table  1.  Averaged  form  drag,  velocity,  power,  and  bulk  drag  coefficient  [(4)].  Estimates  of  Pform  were  limited  to  low-frequency 
motions  (Fig.  11).  The  magnitudes  of  Ki  and  Oi  tidal  velocities  estimated  by  Teague  et  al.  (2013)  are  also  included  in  the  table.  The  values 
inside  the  brackets  are  95%  confidence  limits,  and  the  values  inside  the  curly  brackets  are  minimum  and  maximum  of  the  record.  The  term 
Cd  for  the  diurnal  tides  was  computed  from  velocity  estimates  given  in  Teague  et  al.  (2013).  The  estimated  bulk  drag  coefficient  for 
semidiurnal  tides  by  Warner  et  al.  (2013)  is  given  in  the  bottom  row. 


Variable 

Diurnal  tides  (Ki  and  Oi) 

Subtidal:  4-17-day  variability 

Dform  (Nm-1)  (Fig.  9) 

1500  [1400, 1600]  {350,  4070} 

2326  [1990,  2660]  {352,  4550} 

Uu(cms 

7-11.4  [may  contain  inertial  motions] 

4.5-8.4 

Tidal  analysis 

3-5 

— 

(Teague  et  al.  2013) 

Pform  (W m  );  Pform  ^  Uu 

— 

-85  [-93,  -76] 

(upper-layer  velocity,  P[/  at  Ml) 

Pform  X  U  (depth-averaged 

— 

-73  [-77,  -69] 

velocity;  U  at  Ml) 

Bulk  drag  coefficient  [(4)]:  Cu 

180  X  10“^  {(150-220)  X  10"^} 

36  X  10^^  [(21,  51)  X  10^^]  {(5-180)  X  10^^} 

Cu  from  Fig.  12 

— 

26  X  10“^ 

Cd  (Warner  et  al.  2013) 

Semidiurnal  tides:  90  X  10“^ 

— 

form  drag  can  be  a  linear  function  of  flow  speed  as  found 
by  the  linear  theory  of  lee  waves  generated  by  bottom 
topography  (e.g.,  Baines  1995).  A  formation  of  hydro¬ 
static  lee  waves  and  the  associated  form  drag  for  the 
nondimensional  height  of  the  bank,  /z„  <  1,  and  the 
nondimensional  width,  >  1,  have  been  studied  ex¬ 
tensively  (e.g.,  Baines  1995).  Queney  (1948)  calculated 
the  wave  drag  (Zlwave)  of  a  mountain  terrain  profile 
(“Witch  of  Agnesi”)  with  constant  stratification  and 
constant  background  velocity  provided  /z„  <  1.  In  the 
following,  we  computed  Zlwave  as  a  function  of  time  us¬ 
ing  the  observed  velocity  Uu  and  the  background 
stratification  A,  where 


Fig.  11.  Time  series  of  (a)  upper  50-m  layer  velocity  t/^/  at  MI¬ 
MS,  (b)  Dform  for  the  subtidal  band,  and  (c)  Pform  at  M1-M5.  Thick 
black  lines  in  (a)  and  (c)  are  Uu  and  Pform  at  Ml. 


(16) 

The  minus  sign  in  (16)  is  consistent  with  the  definition  of 
the  form  drag  given  in  (1).  The  term  N  was  estimated 
between  depths  of  15  and  65  m  from  density  time  series 
at  S2,  S3,  and  S4.  The  value  was  taken  as  the  spatially 
averaged  height  of  the  bank  above  the  60-m  bathymetry 
contour  and  is  about  46  m  (±5-m  standard  deviation). 
We  could  not  calculate  N  at  Ml,  so  N  at  S2  was  used.  The 
bin-averaged  P>wave  is  plotted  as  a  function  of  Zlform  in 


Fig.  12.  Form  drag  per  unit  area  Pform/^o  for  the  subtidal  band 
plotted  against  Fb{=—^-^Pq\Uu\Uu)  for  velocities  at  M1-M5. 
The  thick  solid  line  denotes  least  squares  fit  for  Fu  ^  -10.  The 
slope  of  the  line  represents  the  bulk  drag  coefficient  (4)  and  is 
about  26  X  10“^.  The  error  bars  (95%  confidence  levels)  are 
plotted  only  for  Ml. 
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Fig.  13.  Zlform  for  subtidal  motions  vs  the  estimated  wave  drag 
from  the  linear  theory  [(16)]  at  Ml  along  with  95%  confidence 
levels  (vertical  lines). 


Fig.  13.  Both  Zlwave  and  F>form  have  similar  directions, 
but  the  averaged  magnitudes  of  F>wave  are  smaller  than 
^form-  This  could  be  influenced  by  the  terrain-shaped 
Witch  of  Agnesi  used  in  (16),  which  is  different  from  the 
three-dimensional  pear-shaped  EFGB.  There  was  no  ap¬ 
parent  relationship  between  wave  and  form  drags,  espe¬ 
cially  for  negative  T>form  (Fig-  13).  The  correlation 
between  Dwave  and  T>form  (Fig.  13)  is  less  convincing 
compared  to  the  quadratic  drag  approximation  shown  in 
Fig.  12.  The  observed  parameter  range  reflects  a  block¬ 
ing  flow  {hn  ^1),  which  is  significantly  different  from 
Queney’s  nonblocked  flow  solution  (/z„  <  1).  We  discuss 
the  applicability  of  linear  theory  in  the  discussion  section. 

5.  Estimates  of  frictional  drag 

We  compute  the  frictional  drag  coefficient  Cj  [(13)] 
and  the  roughness  length  zo  [(15)]  over  the  bank  by 
combining  measured  TKE  dissipation  rates  and  currents 
near  the  bottom.  The  spacing  between  VMP  profiles  is 
approximately  200  m  on  average,  much  larger  than  in¬ 
dividual  coral  features.  Hence,  our  measurements  can 
be  best  utilized  to  estimate  the  range  of  roughness  over 
the  top  of  the  bank  rather  than  a  precise  value  at  an 
individual  location.  The  measured  values  of  velocity  and 
TKE  dissipation  rate  s  were  interpolated  to  a  grid  with 
0.5-m  vertical  spacing.  The  bottom  boundary  layer 
depth  where  density  changed  by  0.01  kg  m~^  from  the 
near-bottom  density  was  about  4  m  from  the  bottom  for 
nearly  75%  of  the  VMP  profiles  collected  over  the  bank. 
The  inertial  boundary  layer  is  typically  less  than  20%  of 


the  boundary  layer  depth  and  is  expected  to  be  less  than 
1  m.  Dissipation  measurements  were  averaged  over  the 
first  4  m  from  the  bottom  (or  approximately  over  the 
bottom  boundary  layer),  and  five  sequential  profiles  were 
combined,  resulting  in  spacing  on  the  order  of  1  km.  The 
maximum  likelihood  estimator  method  (Baker  and 
Gibson  1987)  was  used  to  calculate  the  expected  value 
and  95th  percentile  confidence  limits.  Velocities  were 
also  averaged  in  time  using  five  adjacent  profiles,  and  any 
measurements  more  than  7.5  m  above  the  bottom  were 
discarded.  The  frictional  velocity  [(14)]  was  calculated 
from  s  as  described  above,  and  errors  are  propagated  to 
the  new  term.  The  near-bottom  mean  velocity  U  was 
calculated  in  two  different  ways:  (i)  using  the  smoothed 
velocity  at  the  deepest  available  depth  (do),  and  (ii)  by  a 
linear  extrapolation  of  velocities  to  3  m  above  the  bottom 
using  the  deepest  two  bins.  These  two  methods  are  re¬ 
ferred  to  as  l7i  and  U2,  respectively.  The  two  techniques 
provide  a  better  sense  of  the  range  of  U,  as  the  ship 
AD  CP  typically  did  not  measure  velocities  within  the 
inertial  layer,  with  the  deepest  bin  typically  5-7  m  above 
the  bottom.  Both  versions  of  U  were  applied  to  calculate 
values  of  zo,  with  errors  in  zo  propagated  from  and  U. 
Calculations  of  Cd  [(13)]  are  based  on  a  reference  height 
at  Za  =  3  m  above  the  bottom,  which  lies  within  the 
bottom  boundary  layer.  The  minimum  for  zo  was  set  to 
10“^  cm,  as  this  is  the  value  for  a  muddy  bottom  (Deacon 
1953),  and  lower  values  are  physically  unrealistic. 

Habitat  composition  maps  show  algal  nodules  cover¬ 
ing  most  of  the  bank,  with  the  coral  reef  zone  (55% 
coverage)  over  the  shallower  regions  (Fig.  1).  The  ex¬ 
pectation  is  that  the  bank  will  contain  a  range  of  values, 
including  both  smooth  areas  and  rough  (coral  coverage) 
regions.  Prior  estimates  indicate  that  the  roughness 
lengths  over  corals  are  expected  to  be  several  centime¬ 
ters  (Rosman  and  Hench  2011),  much  larger  than  what 
was  seen  at  the  moorings  located  just  off  and  on  top  of 
the  bank.  Calculated  roughness  lengths  [(15)],  measured 
dissipation  rates,  and  bottom  velocities  along  the  two 
transects  shown  in  Fig.  14  contain  significant  variability, 
partly  due  to  limited  sampling  in  both  space  and  time. 
Dissipation,  velocity  magnitude,  and  frictional  velocity 
(multiplied  by  10)  are  shown  in  the  first  two  panels  of 
Fig.  14  for  T1  and  T2  (the  transects  highlighted  in 
Figs.  2a,c).  The  third  panel  shows  two  estimates  of  the 
roughness  length  zo,  with  95th  percentile  confidence 
levels  calculated  from  the  near-bottom  velocity  U 
alone  (black  line)  and  with  the  mean  and  standard 
deviation  of  zo  derived  using  Ui  and  U2  (red  line). 
Along  both  T1  and  T2,  which  pass  over  the  shallow 
coral  zone,  the  roughness  lengths  have  values  between 
1  and  10  cm  (Fig.  14),  as  expected  over  corals  (Rosman 
and  Hench  2011). 
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Fig.  14.  Near-bottom  TKE  dissipation  rate  e,  near-bottom  velocity  U,  friction  velocity  (XlO),  roughness  length 
Zo,  and  bathymetry  at  transects  (a)  T1  and  (b)  T2.  The  shaded  regions  in  the  first  three  panels  are  the  95th  percentile 
errors  for  each  parameter.  For  zq,  both  the  value  calculated  using  only  the  near-bottom  velocity  (black  line)  and  the 
mean  and  standard  deviation  of  the  three  Zq  values  obtained  using  Ui  and  U2  (red  line  with  error  bars)  are  shown.  The 
horizontal  line  is  at  1  cm,  the  lowest  reported  value  for  coral  roughness  in  prior  studies  (Rosman  and  Hench  2011). 


The  calculated  roughness  lengths  [(15)]  and  the  as¬ 
sociated  drag  coefficients  [(13)]  over  the  EFGB  vary 
spatially  (Fig.  15).  Given  the  transect  range  of  the  ship 
AD  CP,  estimates  of  zo  and  Cd  are  exclusively  over  the 
top  and  sides  of  the  bank.  To  get  values  over  the  sandy 
bottom  next  to  the  bank,  estimates  of  zo  at  the  long-term 
mooring  locations  were  made  using  times  when  the 
dissipation  transects  came  within  approximately  500  m 
of  a  mooring.  Measured  U  was  1-5  m  off  the  bottom, 
depending  on  the  mooring  configuration.  M4,  M5,  and 
M6  had  close  transects  and  an  AD  CP  bin  within  5  m  of 
the  bottom.  M5  is  near  the  top  the  bank,  but  was  placed 
in  an  area  free  of  coral  in  order  to  avoid  damage  to  the 
reef.  Estimates  for  zo  are  0.02  cm  at  two  mooring  loca¬ 
tions  (M4  and  M6)  and  0.2  cm  at  the  third  (M5).  Some 
uncertainty  is  expected  because  we  do  not  have  dissi¬ 
pation  time  series  at  these  locations,  but  the  values  are 
consistent  with  other  estimates  over  a  sandy  bottom 
(Reidenbach  et  al.  2006).  Distribution  functions  of  both 
Zo  and  Cd  deviate  substantially  from  the  Gaussian  dis¬ 
tribution  and  are  close  to  the  lognormal  distribution. 
The  values  of  zo  based  on  243  estimates  varied  between 


0.001  and  68  cm  with  an  arithmetic  mean  of  3.62  cm 
(Figs.  14, 15;  Table  2).  The  values  of  Cd  varied  between 
1  X  10”^  and  75  X  10”^,  with  an  arithmetic  mean  of  6.2  X 
10“^  (Fig.  15;  Table  2). 

The  range  of  frictional  regimes  over  the  bank  can  be 
seen  in  Fig.  16,  which  shows  Cd  and  the  roughness  Rey¬ 
nolds  number  Re^:,  defined  as  Re^:  =  zoujv.  For  pipe 
flows,  Re^  shows  hydraulically  smooth  regimes  for 
Re^  <  5  and  completely  rough  regimes  for  Re^i-  >70 
(e.g.,  Schlichting  1962).  A  drag  coefficient  of  2  X  10“^  is 
indicated  by  a  black  line,  as  this  is  the  approximate  value 
at  the  mooring  locations  and  represents  a  hydraulically 
smooth  bottom.  The  measurements  comprise  both  values 
typical  for  a  smooth  bottom  as  well  as  values  similar  to 
those  measured  previously  over  corals  (i.e.,  values  of  9- 
150  X  10“^  found  by  Reidenbach  et  al.  2006). 

6.  Discussion 

a.  Form  drag,  power,  and  bulk  drag  coefficient 

Analysis  of  a  6-month-long  time  series  of  pressure 
over  the  EFGB  revealed  that  Dform  was  generated  by 
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Fig.  15.  All  values  of  (a)  zo  (cm)  and  logio  values  of  (b)  bottom  drag  coefficient  (Q)  over  the  bank. 


multiple  time-scale  processes.  The  magnitude  of  T>form 
at  the  diurnal  tidal  band  is  about  1500  Nm“^  for  tidal 
velocities  of  about  3-5  cm  s“^  (Table  1),  and  the  inertial 
drag  was  a  significant  component  of  the  total  form  drag 
for  this  band.  Apart  from  T)form  at  tidal  oscillations,  we 
found  large  form  drags  at  subtidal  bands  with  time  scales 
of  4-6-,  10-11-,  and  16-17-day  periods  (Figs.  9,  10)  for 
depth-averaged  currents  less  than  lOcms”^.  For  these 
low-frequency  bands,  the  inertial  drag  was  negligible. 
Values  of  T>form  for  diurnal  tides  and  4-17-day  oscilla¬ 
tions  were  similar.  Form  drags  resulting  from  semidiurnal 
tides  are  well  documented  in  previous  studies  (Nash  and 
Mourn  2001;  Warner  et  al.  2013),  but  there  are  no  quan¬ 
titative  estimates  of  the  form  drag  for  subtidal  motions 
such  as  4-17-day  oscillations  observed  in  the  GoM.  Our 
analysis  demonstrates  that  the  form  drag  resulting  from 
low-frequency  currents  is  an  important  flow  retardation 
mechanism  even  in  the  presence  of  the  large  frictional 
drag  associated  with  coral  reefs. 

The  6-month  averaged  power  removed  from  the  low- 
frequency  motions  by  the  form  drag  was  85  W  m“^  (Table 
1).  Jarosz  et  al.  (2014,  manuscript  submitted  to  /.  Geophys. 
Res.)  computed  the  along-transect  and  depth-integrated 
turbulent  kinetic  energy  dissipation  rate  E  over  the  bank 
(transects  shown  in  Fig.  2a),  where  E  =  dx.  The 

quantity  E  varied  between  6.5  and  40Wm“^,  and  the 
space-time-averaged  E  over  seven  transects  was  about 
23  Wm“^.  The  short-term  estimates  of  TKE  dissipation 
are  a  factor  of  2-A  smaller  than  the  6-month  averaged 


Fform-  These  differences  are  not  surprising,  since  winds 
and  background  currents  were  weak  during  the  lOP.  The 
6-month  averaged  bulk  drag  coefficient  Cd  over  the  bank 
can  vary  between  0.03  and  0.22,  which  is  3-35  times  larger 
than  the  estimated  frictional  drag  coefficient  Q,  apart 
from  the  fact  that  the  frictional  drag  coefficient  over  the 
bank  is  also  large  due  to  the  coral  habitats  (Tables  1,  2). 

h.  Parameter  dependence  and  linear  theory 

We  computed  the  nondimensional  parameters  hn,  an, 
and  K  (defined  in  section  1)  by  replacing  Nq  by  N  and  Uq 
by  Ujj-  We  noted  for  flow  over  the  EFGB  that  hn  (Fig.  17) 
and  K  were  greater  than  1  for  most  occasions,  except  for 
when  the  stratification  became  weak  and  the  mixed  layer 
was  deep,  as  observed  during  late  fall  frontal  passages. 
Here  N  is  estimated  from  the  6-h  low-pass  filtered  po¬ 
tential  density  at  15  and  65  m  at  S2,  and  f/j/  is  the  6-h  low- 
pass  filtered,  upper  50-m,  depth-averaged  flow  speed  at 

Table  2.  Bottom  roughness  zo  and  frictional  drag  coefficient  Q 
estimated  over  the  EFGB.  Maximum,  minimum,  arithmetic  mean, 
and  median  were  estimated  from  243  observations.  The  maximum 
likelihood  estimator  (MLE)  was  computed  by  assuming  lognormal 
distributions  for  zo  and  C^,  where  MLE  =  exp(m  +  s^/2);  m  and  5 
are  the  arithmetic  mean  and  standard  deviation  of  /n(variable) 
(e.g..  Baker  and  Gibson  1987). 


Variable 

Minimum 

Maximum 

Mean 

Median 

MLE 

Zo  (cm) 

0.001 

68 

3.6 

0.24 

22 

Cd 

0.001 

0.075 

0.0062 

0.0032 

0.0056 
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Fig.  16.  The  term  Cj  as  a  function  of  frictional  Reynolds  number 
Re^.  The  horizontal  line  is  at  0.002,  the  approximate  value  of  Cj  at 
the  mooring  locations. 


Ml.  For  /zo  =  46  m,  D  =  100 m,  L  =  6  km,  and  =  1  km, 
holL  is  0.017  and  hJD  is  0.46.  The  nondimensional  pa¬ 
rameters  hn,  K,  and  calculated  from  the  time-averaged 
velocity  Uu  and  the  buoyancy  frequency  N  are  about 
5.6  (±3.5),  3.8  (±2.6),  and  125  (±78),  respectively,  where 
the  time-averaged  Uu  and  N  are  0.15  (±0.07)  ms“^,  and 
0.023  (±0.005)  s“^,  respectively.  The  standard  deviations 
are  given  in  parentheses.  The  height  parameter  Hq/D 
shows  that  the  bank  occupies  nearly  50%  of  the  water 
column,  and  values  of  nondimensional  parameters  (tz„  > 
1,  /z„  >  1,  >  1,  and  DIL  <C  1)  suggest  highly  nonlinear, 

hydrostatic,  high  drag  state  flow  with  upstream  blocking 
(e.g.,  Baines  1979;  Hunt  and  Snyder  1980;  Vosper  et  al. 
1999).  The  relation  >  1  indicates  that  the  lateral  scale 
of  the  bank  is  much  larger  than  the  distance  traveled  by 
the  flow  during  a  buoyancy  period. 

Garrett  and  Kunze  (2007)  have  suggested  extra  pa¬ 
rameters  including  a  steepness  parameter,  = 
/i:/zo[(ct)^  — /^)/(N^  —  and  a  tidal  excursion 

kUu/(Ji>  for  oscillatory  flows  over  the  bank  with  fre¬ 
quency  (o  (>/)  and  horizontal  wavenumber  k.  For  flows 
over  the  EFGB,  the  steepness  parameter  exists  for  the 
Ki  tide  with  coKi  (=7.29  X  10“^  s“^)  >  /  (=6.84  X 
10~^s~^),e„  =  10"^and/:f/[//mKi  =  2.0for/:  =  1/ao  =  1  X 
10“^  m“^.  The  relation  >  1  suggests  that  the  topo¬ 
graphic  slope  is  supercritical  with  respect  to  the  tidal  ray 
slope,  and  kUuloj  =  2  suggests  that  quasi-steady  lee 
waves  can  be  generated  in  the  region. 

Formation  of  hydrostatic  lee  waves  and  the  associated 
form  drag  for  >  1  and  /z„  <  1  have  been  studied 


20 


18 


200  250  300 

Year  Day 


Fig.  17.  Nondimensional  bank  height  hn  at  Ml.  Fluctuations  less 
than  10  h  were  filtered  out  before  plotting.  The  dashed  line  rep¬ 
resents  hn  =  l. 

extensively  (e.g.,  Baines  1995).  As  discussed  in  section  4e, 
Dwave  and  T>form  havc  similar  directions,  but  the  magni¬ 
tude  of  T>wave  is  Smaller  than  the  measured  form  drag 
(Fig.  13).  The  dynamics  of  the  blocking  flow  (hn  >  1) 
observed  here  was  significantly  different  from  Queney’s 
(1948)  nonblocked  flow  solution.  A  theoretical  de¬ 
scription  of  the  flow  for  >  1  and  /z„  >  1  is  not  well 
established.  If  the  wave  drag  [(16)]  can  be  approximated 
as  the  form  drag  defined  in  (3)  and  (4),  then  drag  co¬ 
efficients  are  functions  of  nondimensional  parameters 
such  as  hn  and  ho/lo,  that  is,  Cdf  =  7r/z„/4  and  Co  = 
7Thnl4(hollo).  Vosper  et  al.  (1999)  show  that  the  Cdf  in¬ 
creases  by  a  factor  3  as  hn  increases  from  1  to  4  and  is 
consistent  with  Cdf  ^  TrhJ^.  This  linear  relationship 
between  Cdf  and  hn  was  explored,  but  the  error  bars 
were  large,  and  no  functional  form  was  found. 

Jarosz  et  al.  (2014,  manuscript  submitted  to/.  Geophys. 
Res.)  reported  upstream  stagnation  flow,  flow  around  the 
bank,  vortex  formation,  and  bottom  flow  reversals  on  the 
lee  side  of  the  EFGB.  They  noted  that  their  findings  were 
qualitatively  similar  to  the  laboratory  experiments  (e.g., 
Baines  1979;  Hunt  and  Snyder  1980;  Vosper  et  al.  1999). 
The  dynamical  parameters  (a„  >  1,  /z„  >  1,  >  1,  and 

DIL  <C  1)  and  the  observations  of  Jarosz  et  al.  (2014, 
manuscript  submitted  to  /.  Geophys.  Res.)  suggest  that 
at  least  part  of  the  measured  form  drag  was  caused  by 
pressure  anomalies  associated  with  flow  separation  and 
vortices  generated  on  the  lee  side.  The  strongest  flow 
over  the  bank  was  associated  with  near-inertial  waves 
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(Teague  et  al.  2013);  hence,  flow  separation  and  eddy 
shedding  could  occur  within  the  inertial  cycle  on  the  lee 
side  of  the  bank  (especially  over  the  steep  slopes),  apart 
from  the  low-frequency  motions.  The  Rossby  number  Uu! 
fL  is  about  0.35  (±0.14)  for  Uu  =  0.15ms“\/=  6.84  X 
10“^  s“^,  and  L  =  6  km,  and  therefore  rotation  can  also 
be  a  factor  in  the  local  dynamics  (e.g..  Smith  1979). 

c.  Frictional  drag  coefficients  associated  with  corals 

The  frictional  drag  coefflcient  over  the  bank  is  large 
due  to  high  roughness  elements  associated  with  coral 
reefs.  The  observations,  which  were  utilized  to  estimate 
this  coefflcient,  were  made  when  winds  were  weak  and 
surface  waves  were  small.  Consequently,  surface  wave- 
induced  bottom  currents  could  not  influence  the  bottom 
TKE  dissipation  rate.  The  estimated  roughness  lengths 
reflect  the  true  roughness  lengths  of  small-scale  bathy¬ 
metric  features  over  the  EFGB.  The  coral  coverage  is 
not  continuous  over  the  bank,  and  thus  zo  had  a  range  of 
values  from  those  representing  a  smooth  bottom  with 
a  roughness  length  of  0.001  cm  to  those  for  relatively 
rough  regimes  with  zo  as  high  as  68  cm.  Nearly  Vs  of  our 
estimates  were  greater  than  1  cm,  and  the  average  drag 
was  about  6  X  10”^.  These  higher  roughness  lengths 
correspond  to  high  drag  coefficients  and  large  frictional 
Reynolds  numbers.  The  impact  of  high  bottom  frictional 
forces  can  be  important  because  the  strength  of  the  bot¬ 
tom  boundary  layer  stress  can  change  lee-wave  forma¬ 
tion,  decrease  the  form  drag,  and  reduce  the  overall 
momentum  loss  to  bottom  obstacles/banks  compared  to  a 
typically  assumed  free-slip  boundary  (Skyllingstad  and 
Wijesekera  2004).  Pratt  (1986)  reported  that  hydraulic 
controlled  flows  over  straits  and  sills  with  friction  can 
alter  the  dynamics  of  the  flow  such  as  shifting  of  control 
points  from  a  sill  to  downstream  locations.  Although 
these  modeling  studies  do  not  mimic  the  EFGB  observa¬ 
tions,  the  observed  large  bottom  stress  over  the  bank  is  a 
factor  to  consider  in  the  parameterization  of  the  form 
drag. 

7.  Summary  and  conclusions 

Hydrographic,  velocity,  microstructure,  and  bottom 
pressure  measurements  were  used  to  quantify  the  mag¬ 
nitude  and  temporal  variability  of  the  form  drag,  rough¬ 
ness  lengths,  and  associated  frictional  drag  coefficients 
over  the  EFGB  bank.  The  EFGB  is  a  rough  topographic 
feature  located  on  the  Louisiana-Texas  shelf  approxi¬ 
mately  190  km  southeast  of  Galveston,  Texas.  The 
EFGB,  about  6  km  wide  and  10  km  long,  is  located  at  the 
shelf  edge  in  100  m  of  water  depth  with  a  peak  rising  to 
about  18  m  below  the  sea  surface.  The  area  encompassed 
by  the  100-m  bathymetric  contour  of  the  bank  is  about 


40  km^.  The  bank  is  one  of  the  northernmost  tropical 
coral  reefs.  The  bank  includes  both  coral  coverage  and 
sandy  bottoms. 

Nearly  6-month-long  bottom  pressure  records  (June- 
December  2011)  showed  significant  semidiurnal  (M2)  tide, 
diurnal  (Ki  and  Oi)  tides,  and  4-17-day  variability  over 
and  around  the  bank.  Note  that  tidal  analysis  conducted 
by  Teague  et  al.  (2013)  showed  that  the  barotropic  tidal 
currents  were  small  (Ki,  Oi  speeds  <  5  cm  s“^)  compared 
to  inertial  currents  that  were  over  15  cm  s~^.  The  M2  tidal 
currents  were  less  than  2  cm  s“^,  and  the  depth- averaged 
subtidal  currents  were  about  4-9cms“^. 

The  estimated  form  drag  from  bottom  pressure  anom¬ 
alies  across  the  bank  showed  variability  on  multiple  time 
scales.  Here  T>form  was  estimated  from  three  pressure 
measurements  over  and  around  the  bank  while  assuming 
that  the  local  horizontal  pressure  gradient  was  a  constant 
between  two  observational  sites  (e.g.,  Hafner  and  Smith 
1985).  It  is  likely  that  T>form  was  underestimated  by 
a  factor  of  2.  Therefore,  our  estimate  of  T>form  could  be 
treated  as  a  lower  bound  of  the  expected  value.  The  drag 
resulting  from  diurnal  tidal  motions  and  4-17-day  oscil¬ 
lations  had  magnitudes  of  about  2000  N  m“^.  Coefficients 
of  the  time-averaged  bulk  drag  Cd  ranged  from  0.03  to 
0.22,  which  was  a  factor  of  5-35  larger  than  the  averaged 
frictional  drag  coefficient  (Tables  1,  2).  The  term  Pform 
was  estimated  as  a  product  of  the  velocity  and  the  form 
drag  for  low-frequency  motions.  The  enumerated  power 
loss  per  unit  length  at  subtidal  motions  was  about 
85  W  m“^.  The  total  power  loss  over  the  bank,  PformAo/Zo, 
was  —0.37  MW,  where  Aq  ^  40  km^  and  Iq  =  9  km. 

Physical  parameters  describing  the  dynamics  of  flows 
over  the  EFGB  have  a  range  of  values  representing  a 
highly  nonlinear,  hydrostatic  flow  regime  with  upstream 
blocking,  where  /z„  >  1,  >  1,  and  an  —  0(100).  For 

these  parameters,  there  is  no  theoretical  prediction  of  the 
form  drag.  The  wave  drag  based  on  linear  theory  of  the 
flow  over  a  bank  with  hn  and  the  measured  drag  for 
subtidal  flows  had  similar  directions,  but  the  linear  theory 
underpredicted  the  observations.  The  quadratic  drag  law 
appears  to  have  a  better  agreement  with  observations. 

Using  measurements  of  the  velocity  and  estimated 
TKE  dissipation  rates  along  multiple  transects  crossing 
the  bank,  the  range  of  roughness  lengths  zo  and  bottom 
frictional  drag  coefficients  were  computed  (Table  2). 
These  measurements  of  the  velocity  and  dissipation  over 
the  bank  provide  the  first  estimates  of  frictional  drag  co¬ 
efficients  and  roughness  lengths  over  the  EFGB.  Over  the 
top  of  the  bank  the  bottom  roughness  lengths  encompass 
those  expected  for  a  smooth  bottom  with  zo  of  0.001  cm  as 
well  as  higher  roughness  values  over  corals.  Given  that 
the  coral  coverage  is  not  continuous,  a  range  of  values 
representing  smooth  to  relatively  rough  regimes  was 
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expected.  The  range  of  zo  was  from  0.001  to  68  cm,  with 
mean  and  median  values  of  3.6  and  0.24  cm,  respectively 
(Table  2).  Nearly  Vs  of  the  roughness  estimates  were 
greater  than  1  cm.  The  mean  and  median  of  the  frictional 
drag  coefficient  were  about  6  X  10”^  and  3  X  10”^, 
respectively  (Table  2).  These  findings  are  comparable 
with  previous  measurements  over  other  coral  reefs  (e.g., 
Reidenbach  et  al.  2006). 

Large  roughness  elements  can  impact  the  hydraulic 
flow  control  as  well  as  coral  reef  ecology.  Locations  with 
coral  coverage  have  increased  drag,  leading  to  more 
turbulence  and  mixing  at  these  locations  relative  to  re¬ 
gions  without  corals.  This  turbulence  in  turn  increases 
the  delivery  of  nutrients  and  particulates  to  the  corals. 
Additionally,  it  helps  to  disperse  larvae  during  spawning 
events,  generating  a  positive  feedback  cycle  that  sustains 
the  health  of  the  coral  community.  The  magnitude  of  the 
bottom  stress  over  the  bank  can  modify  the  separation  of 
flow,  lee-wave  structure  and  strength  of  the  form  drag; 
hence,  it  is  a  factor  to  consider  in  the  parameterization 
of  the  form  drag  in  numerical  models.  Consequently,  our 
observations  and  analysis  of  the  hydrographic,  velocity, 
microstructure,  and  pressure  measurements  suggest  that 
the  EFGB  is  a  “hot  spot”  of  mixing  on  the  shelf  of  GoM. 
The  analysis  further  demonstrates  that  the  form  drag  re¬ 
sulting  from  low-frequency  currents  over  an  isolated  bank 
on  the  continental  shelf  is  an  important  physical  process 
that  must  be  parameterized  to  represent  a  wide  range  of 
flow  states. 
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APPENDIX 

Estimation  of  the  Form  Drag  Based  on  Numerical 
Simulations:  Sensitivity  Study 

The  form  drag  over  an  underwater  topographic  fea¬ 
ture  is  estimated  from  bottom  pressure  fluctuations  and  a 
horizontal  slope  of  the  bathymetry,  as  indicated  in  (1). 


The  accuracy  of  the  estimated  drag  depends  on  the 
number  of  pressure  measurements  and  the  location  of 
the  pressure  measurements.  A  sensitivity  study  was 
conducted  to  examine  the  minimum  number  of  pres¬ 
sure  measurements  required  to  calculate  a  reasonable 
estimate  of  the  form  drag.  Here  we  used  numerically 
simulated  pressure  fields  over  a  sloping  headland,  TTP 
in  Puget  Sound,  Washington,  developed  by  Warner 
(2012)  as  part  of  her  Ph.D.  dissertation.  The  TTP 
model  was  based  on  the  Regional  Ocean  Modeling 
System  (ROMS),  and  details  of  the  model  setup  is 
given  by  Warner  (2012).  The  study  helps  us  to  quantify 
some  uncertainties  in  an  estimation  of  the  form  drag 
over  the  EFGB  based  on  three  pressure  measure¬ 
ments.  The  ROMS  results  are  an  approximation  to  the 
real  observations  due  to  numerical  limitations  given  by 
Warner  (2012).  It  is  also  noted  that  the  shape  and  size 
of  TTP  are  different  from  the  EFGB  and  that  the  lo¬ 
cations  of  Ppods  during  MORT  might  have  missed 
large  pressure  fluctuations.  Therefore,  the  results  dis¬ 
cussed  below  have  more  uncertainties  than  the  model 
results  alone. 

A  time  series  of  modeled  bottom  pressure  anoma¬ 
lies  along  the  TTP  ridge  was  used  to  compute  the 
form  drag  based  on  (1).  The  spatial  resolution  of  the 
model  was  50.5  m.  The  ridge  is  about  2  km  wide  and 
is  nearly  symmetrical  (Fig.  Al).  Gradients  of  the  ridge, 
perpendicular  to  the  ridge  axis,  are  strongest  on  the 
midsections  of  the  ridge  with  the  maximum  on  the 
right  side  and  the  minimum  on  the  left  side.  Total 
pressure  anomalies  [(8)]  including  external,  internal, 
and  inertial  components  at  three  different  locations 
(near  bottom  and  top  of  the  ridge)  are  shown  in 
Fig.  A2. 

We  computed  the  form  drag  using  the  modeled 
pressure  along  the  ridge  and  the  measured  slope  with 
50.5-m  spatial  resolution  and  referred  to  it  as  the  com¬ 
plete  solution  Dq. 

^0  =  “  L  (Pb),  and  \x  =  50.5  m,  (Al) 

where  I  is  the  total  number  of  pressure  measurements 
used  in  the  integration,  and  X3  —  Xi  =  I  Ax  =  2.1km 
(Fig.  Al).  We  computed  Dform  by  approximating  the 
pressure  field  based  on  three  pressure  measurements  at 
Xi,  X2,  and  X3,  while  mimicking  the  observational  setup 
over  the  EFGB  during  the  MORT. 

a.  Method  I 

The  first  estimate  was  made  by  using  mean  values  of 
the  pressure  and  slope  over  the  bank,  and  the  resulting 
drag  is 
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Fig.  A1.  (a)  Cross  section  of  the  ridge  transect  based  on  the  model  bathymetry.  Spatial 
resolution  is  50.5  m.  The  crosses  are  the  data  points  used  in  the  analysis.  The  bullets  represent 
locations  of  pressure  time  series  plotted  in  Fig.  A2.  (b)  Horizontal  gradient  of  ridge  topography, 
dhidx  across  the  ridge  (from  Warner  2012). 


^1  ^ml 


(A2) 


where  Pmi  =  +  pb(X2)]/2  and  Pmi  =  [Pb(X2)  + 

Pb(X2)]/2.  The  angle  brackets  are  the  spatial  means  over 
the  slope. 

b.  Method  II 

The  second  approximation  assumed  that  the  pressure 
distribution  was  linear  over  the  slopes.  We  interpolated 
pressure  linearly  between  Xi  and  X2  representing  the 
left  side  of  the  bank  and  pressures  between  X2  and  X3 
representing  the  right  side  of  the  bank.  The  resulting 
drag  based  on  the  linearly  interpolated  bottom  pressure 
Pi  with  fully  resolved  bathymetric  slope  {dhidx)  is 

f  (/>,), (A3) 


c.  Method  III 


Here  the  form  drag  was  approximated  by  using 
a  mean  value  of  the  pressure  and  the  maximum  slope 
over  the  bank,  and  the  resulting  drag  is 


^3  ^ml 


maxi 


[X^-XJ  +  F, 


m2 


[X3 


max2 


X2], 

(A4) 


The  estimated  drag  over  the  ridge  from  the  three 
different  methods  is  illustrated  in  Fig.  A3.  The  drag  re¬ 
sulting  from  all  three  methods  have  the  same  sign,  even 
though  the  magnitudes  of  the  complete  solution  {Dq) 
and  method  III  (D^)  are  larger  than  Di  [(A2)]  and  D2 
[(A3)].  The  terms  Di  and  D2  have  similar  magnitudes, 
but  their  maximum  magnitudes  are  about  2  times 
smaller  than  the  maximum  of  and  £>3  (Fig.  A3). 
However,  when  averaged  over  the  tidal  cycle,  both  Di 
and  D2  are  about  70%  of  the  complete  solution  Dq. 


where  \  dhldx\^^^-^  and  \dhldx\^^^^  are  the  maximum  slopes 
of  left  and  right  sides  of  the  bank  (Fig.  2b). 


Fig.  A2.  Dynamic  bottom  pressure  fluctuations  at  three  different 
locations:  Xi,  X2,  and  X3.  Locations  of  Xi,  X2,  and  X3  are  marked 
in  bullets  in  Fig.  A1  (from  Warner  2012). 
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Fig.  A3.  Estimated  pressure  drag  from  four  different  methods: 
the  complete  integration  Dq  [(A1)],  mean  pressure  and  mean  lat¬ 
eral  gradient  Di  [(A2)],  linearly  interpolated  solution  D2  [(A3)], 
and  mean  pressure  and  maximum  amplitude  of  lateral  gradient,  D3 
[(A4)]. 


There  are  large  variations  of  the  slope  across  the  ridge, 
and  it  appears  that  the  slope  is  an  important  factor  in  the 
computations. 
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